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ABSTRACT 


Exhaust gas temperature, velocity and composition are 
measured and combustor efficiencies are calculated in a lean 
premixed swirl-stabilized laboratory combustor. These data 
should be valuable in understanding the combustion and fluid 
mechanical processes in such a combustor with an aim to 
reduce pollutant emissions and maintain high combustion 
efficiencies in order to develop new combustor designs for 
practical applications. This work is part of an overall 
combustion research program at Cornell University directed 
towards the characterization of swirling turbulent reacting 
flows . 

The combustor consists of two confined co-axial swirling 
jets. The inner jet is premixed fuel-air mixture while the 
outer jet is air only. The amount of swirl can be varied in 
both the jets. The outer jet can be swirled in the same 
direction (co-swirl) or in the direction opposite (counter- 
swirl) to that of the inner jet. Combustion is stabilized in 
front of a swirl-induced recirculation zone. 

In the first phase of measurements, radial traverses 
measuring temperature, velocity and gas composition are made 
at the exit plane for varying outer swirl conditions using 
propane as fuel. These data are compared with previous 
results for methane firing. In the second phase, 
measurements are made under cold flow (nonreacting) 
conditions to investigate fuel distribution at the inlet and 


exit planes. In addition, data similar to those obtained in 
the first phase are obtained for methane firing for 
comparison with data for propane firing as well as previous 
results . 

Temperature measurements are made with an uncoated 
Pt/Pt-105&Rh fine wire thermocouple. A transverse cylindrical 
probe is used to 'measure velocity as well as to draw gas 
samples. In situ calibration of the probe is found necessary 
for velocity measurements to correct for turbulence effects. 
Samples are analyzed for oxides of nitrogen in a 
chemiluminescence analyzer equipped with an NO^-NO converter. 
Analyses for the major species are done by gas 
chromatography. A novel scheme to significantly reduce the 
analysis time for the gas chromatographic analysis of the 
products of propane combustion is developed. 

Results show significant differences in the radial 
profiles of the data between the co- and the counter-swirl 
cases. Co-swirl cases show evidence of poor turbulent mixing 
across the combustor in comparison to the counter-swirl 
cases. NO^ levels are low in the combustor but substantial 
amounts of CO are present. Combustion efficiencies are low 
and surprisingly constant with varying outer swirl in 
contradiction to previous results under a slightly different 
inner swirl condition. This difference in the efficiency 
trends is expected to be a result of the high sensitivity of 
the combustor to changes in the inner swirl. Combustor 


operation is found to be the same for propane and methane 
fuels . 

A mechanism is proposed to explain the combustor operation 
and a few important characteristics determining combustor 
efficiency are identified. Future experiments to 
substantiate the proposed mechanism are suggestf;d. 
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CHAPTER 1 

I 

INTRODUCTION 

^ Swlrl-Stabllized Combustion 

Swirling flows have been studied for many years because of 
their potential for stabilization of high intensity combus- 
tion processes [1,2]|, Swirl has been used to stabilize sub- 
sonic combustion systems[33 and is proposed for use in super- 
sonic combustion processes as well[4]. The attractive 
features of reacting vortex flows are better flame stability 
resulting from the formation of one or more recirculating 
flow zones, reduced combustion length and increased mixing 
through turbulence generation. These large scale effects 
have been demostrated by numerous experimental studies 
involving both diffusion and premixed flame combustors, 
e .g . I C3> 5> 6, 7, 8]. Unfortunately, the simultaneous presence 
of swirl, shear, turbulence and heat release makes the react- 
ing vortex flow a very complex object of study and as a 
consequence these flows are poorly understood. 

Recirculation in flows with swirl is known to depend on 
the swirl number S, the ratio of the axial flux of angular 
momentum to the axial flux of linear momentum, and the 
Reyolds number of the flowCI]. As the swirl number is 
increased a pronounced deceleration or even a reversal of the 
flow along the flow axis may occur. Below a critical swirl 
number no back flow occurs. Above this critical level of 
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pwirl the phenomenon of vortex breakdown[9 , 10] occurs which 
is characterized by the formation of a free stagnation point 
on the vortex axis, followed by reversed flow in a region of 
limited axial extent. Various forms of vortex breakdown are 
possibleC9 , 1 1 ] . Notable among them are the axisymmetric and 
the spiral forms, both of which appear capable of stabilizing 
combustion[ 1 2] . Continued reasearch is required to under- 
stand and characterize combustion processes in swirling flows 
and to develop combustor designs with such flows. 

One of the present applications of the swirl-combustor is 
in gas turbine engines. It is hardly necessary to stress the 
increasing role of gas turbines in the present day world. 
Gas turbines are the major power systems for military and 
commercial aircraft. Also, gas turbines are used for marine 
propulsion, in the ’Spruance* class of destroyers, for exam- 
ple. A recent application is in electric power generation as 
peaking-power units, mid-range and even as baseload power[133 
units. In the combined cycle concept a gas turbine topping 
cycle generates power and its exhaust is used to provide heat 
to the bottoming cycleCl4]. As a result of the high operat- 
ing temperatures possible with new high temperature alloys 
and ceramic materials, automotive gas turbines are under ser- 
ious consideration[ 1 5] . 

1 . 2 Pollutant Emissions and Control 

With the growth in applications and the widespread use of 
gas turbines, the problem of pollutant emissions from gas 
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turbines also has grown. The need to limit pollutant emis- 
sions from combustion systems is widely recog- 
ni2ed[ 1 6, 17> 18, 19, 203 . The U.S. Environmental Protection 
Agency (EPA) has proposed stringent emission standards for 
new airoraftC21] and stationary gas turbines[22]. The three 
major pollutants from gas turbine combustors are carbon 
monoxide (CO), unburned hydrocarbons (UHC) and oxides of 
nitrogen (NO ) which comprise nitric oxide (NO) and nitrogen 

A 

dioxide (NOg). Soot and sulphur oxides emissions can occur 
under certain conditions but do not currently pose problems 
for gas turbines. 

1.2.1 Unburned Hydrocarbons 

Unburned hydrocarbons (UHC) are a direct indication of 
inefficient combustion. In addition, they are intimately 

connected with oxides of nitrogen in the formation of photo- 
chemical smog in the atmosphere. Higher UHC emissions are 
caused by increasingly lean conditions which lead to lower 
flame temperatures and to thermal quenching of the hydrocar- 
bon breakdown reactions. Thermal quenching can be described 
as the cooling of the gases below the temperature at which a 
given reaction can effectively take place. Sawyer[23] exper- 
imentally determined that the quenching temperature for heavy 
unburned hydrocarbon oxidation is in the range of approxi- 
mately 1500 to 1700°K. Also, large droplets with long burn- 
ing times can contribute directly to increased UHC as they 
may persist into the secondary dilution zone or hit the cool 



combustor liner and remain unreacted. Thus premixed, 
prevaporized combustion in a near stoichiometric mixture is 
likely to reduce or eliminate UHC emissions. 

1.2.2 Carbon Monoxide 

In addition to being highly toxic, carbon monoxide in the 
gas turbine exhaust is also an indication of inefficient com- 
bustion-, It is well accepted that the oxidation of CO to COg 
in combustion systems follows the reaction[2i| ] 

CO + OH = COg + H . (1.1) 

Reaction (1.1), which proceeds rapidly at high temperatures, 
is hindered at low temperatures due to reduced levels of OH. 
At stoichiometric conditions the CO concentrations are kept 
to a minimum in the primary zone of the combustor where high 
flame temperatures are reached and combustion is complete. 
CO increases in rich mixtures due to lack of oxygen and in 
lean mixtures due to thermal quenching by excess air. This 
similarity between HC and CO trends in lean flames is due 
mainly to the tendency of oxidation reactions to increase 
with temperature. Hov/ever, since the CO oxidation is slower 
than hydrocarbon reactions, CO is more prone to thermal 
quenching and persists unreacted for many (5-6) milliseconds 
in fuel lean flames (s^ =0 . 8) [ 1 7] . The quenching temperature 
for CO has been experimentally determined [ 1 73 to be about 
1270 '^K at atmospheric pressure and approximately 1500°K at 3 
atmospheres below which temperature CO can remain unreacted 
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indefinitely if quenching ii fast enough. These 

observations[ 1 73 are consistent with the approximation that 
flame generated trace species such as 0, H, Hg and OH remain 
equilibrated among themselves and decay together as the temp- 
erature falls, but CO is not necessarily equilibrated with 
the others. Since the CO reaction with OH can be slow com- 
pared to the decay of the trace species, CO can remain 
largely unreacted as the others decay. At high temperatures, 
however, CO can occur due to dissociation of COg. Reduction 
in CO emissions can be achieved by uniform, stoichiometric or 
lean premixed combustion with reduced thermal quenching 
rates . 


1.2.3 Oxides of nitrogen 

NO and NOg are highly toxic. In addition, they play an 
important role in smog formation and depletion of ozone in 
the stratosphereC 1 6]. Since, at ground level, an intercon- 
version between NO and NOg is possible in the presence of 
atmospheric oxygen, there has been little interest in spe- 
cific emissions of NO and NOg, and emissions have been mea- 
sured in terms of NO . In the stratosphere the catalytic 

A 

ozone destruction cycle promotes a steady-state balance bet- 
ween the oxides of nitrogen on a relatively short time 
scale[19] regardless of the initial amounts of each. Smog 
formation, however, may well be affected by high NOg levels, 
since ozone does not begin to form until NO has been con- 
verted to NOg in the atmosphere. It has also been suggested 
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that locally high NO^ emissions in the lower stratosphere 
could, under some conditions, result in a more rapid removal 
of via conversion to nitric acid[25l. Thus the relative 
amounts of NO and NOg are important and in recent studies the 
two species are reported separately. 

In combustion systems, NO formation is a prerequisite for 
NOg formation. It is generally accepted that NO formation in 
fuel-lean mixtures follows the Zel'dovich mechanism (thermal 
NO) where fuel-bound nitrogen is insignificant[263 . Forma- 
tion of NO is associated with high temperatures of 1900°K and 
above[23]. Under lean conditions in a combustor NO increases 
with increasing equivalence ratio, following a trend opposite 
that of CO and UHC, due to the increase in temperature with 
increasing equivalence ratio. But even at high temperatures 
NO formation is fairly slow and occurs on a time scale of 
several milliseconds. Hence NO is formed in the primary zone 
where the two main prerequisites for NO formation are met, 
namely, high temperatures and sufficient residence time. 
This understanding suggests two methods of NO control. 
First, the peak temperatures in the combustor should be low- 
ered by methods such as lean premixed combustor operation and 
exhaust gas recirculation. In the case of liquid fuels whose 
individual droplets burn at locally stoichiometric condi- 
tions, prevaporization and premixing of fuel seems necessary. 
Second, the residence time in the primary combustion zone 
should be reduced by the rapid addition of secondary air 
through turbulent mixing. Schefer and Sawyer[26] have pre- 
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sented results for turbulent, lean, premlxed combustion, of 
their own work as well an others which show substantial 
reductions in the levels of emissions over the emis- 
sions from current conventional turbulent diffusion flame gas 
turbine combustors, Gouldin[273 points out that while turbu^ 
lent mixing is desirable both to reduce residence times in 
high temperature zones and to enhance the stabilization of 
the lean combustion suggested above, the intensity and small 
scale of the turbulence required could result in a prohibi- 
tive combustor pressure drop. 

The proceeding discussion applies to NO formed from Ng in 
atmospheric air. Nitrogen is also found in some fuels, espe- 
cially in those derived from coal. Fuel nitrogen conversion 
to NO occurs at lower temperatures than reactions with mole- 
cular nitrogen and may result in larger amounts of NO than in 
pure hydrocarbon flames[28]. It is suggested that flame gen- 
erated nitrogeneous species such as HCN, CN and NH^^ play 
important roles in NO formation in flames of fuels contain- 

A 

ing nitrogen compounds[ 283* Also, in fuel-rich hydrocarbon 
flames, some NO formation proceeds by a mechanism similar to 
that for fuel-bound nitrogen involving the formation of 
intermediate species ’such as HCN and CNL’293. 

The formation of NOg is less completely understood than 
the formation of NO though several explanations have been 
tendered. A consistent mechanism for NOg formation must be 
able to explain the high fractions of NOg that have been mea- 
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sured in or near laminar flame frontsCSO, 31 , 32] as well as in 
turbulent combustorsC26, 33, 3^, 353. The most plausible 
mechanism for NO to NOg conversion invokes Reaction (1.2) in 
quenching zones where the temperature is falling and the free 
radical concentrations exceed their equilibrium 
values[25,31,32,353. 



r 

NO ■ 

f HOg = 

NO 2 + OH 

(1.2) 

FenimoreC 32 ] has 

shown 

that with the assumption 

of local 

equilibrium for 

the 

Reactions 

(1.2), (1.3) and 

(1.4) the 

ratio of NO 2 /NO 

can 

be 

found 

by assuming equilibrium for 


Reaction (1.5) when oxygen atoms are determined by the equi- 
libration of Reaction (1.4). 


°2 * M 

= HO 2 + M 

(1.3) 

H + Og 

= OH + 0 

(1.4) 

NO + 0 

= NO 2 

(1.5) 


Calculations by Anderson, et al.[25] have shown that under 
after burning conditions in a supersonic transport engine, 
super equil ibrium concentrations of HOg can result in unex- 
pectedly large values for NOg/NO through Reaction (1.2). 
This supports Fenimore’s results although in Fenimore's case 
the excess HO^ results from hydrocarbon oxidation rather than 

Im* 

from Reaction (1.3). Thus, Fenimore’s mechanism with excess 
HO 2 due to hydrocarbon oxidation appears to be the most plau- 
sible explanation of NO 2 formation in gas turbines. Further 
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support for this mechanisni ponies from the study by Chen, ct 
al.[36] who showed that when a partially reacted sample con- 
taining unburnt fuel is quenched (either by mixing or in a 
probe), large concentrations of HOg are produced during low 
temperature hydrocarbon oxidation. Reaction (1.2) is there- 
fore favored in quenching zones, and the high NOg/NOj^ frac- 
tions can be attributed to quenching zones in the combustor. 

1.2.^ Pollution Control 

Some conditions under which pollutant emissions can be 
reduced have been mentioned in the preceding subsections. 
Pollution control in gas turbines and other continuous flow 
combustors has been accomplished by trying to achieve these 
conditions primarily by alterations or additions to existing 
designs. Since conditions minimizing NO formation tend to 
maximize CO and UHC emissions, and vice versa, a trade off is 
inevitable. Flame temperatures have been reduced by leaner 
fuel-air mixtures, by steam or water injection into the com- 
bustion zone and by exhaust gas recirculation. Fuel injec- 
tion has been improved and combustor liners have been modi- 
fied to change air flows. LefebvreC 37 3 , in a review of gas 
turbine pollution control, concludes, as do Rudy and ReckCl6] 
and many others, that only advanced and radically new combus- 
tor designs involving concepts such as prevaporized premixed 
combustion, variable geometry and swirl-augmented combustion 
can achieve the required emissions reduction levels. These 
are the very concepts which form the basis for the swirl- 
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stabilized combustion research at Cornell. 

^ • 3 Swirl-Stabilized Combustor Research at Cornell 

Research at Cornell in this trea, going on since the early 
seventies, is a continuing program directed towards the study 
of the interrelationships between chemical kinetics and tur- 
bulent mixing in , a research combustor with premixed pre- 
vaporized fuel and variable swirl. The aim of the overall 
program is to completely characterize the swirling flows in 
the combustor under both reacting and nonreacting conditions 
and perhaps to develop numerical models to predict these 
flows, This characterization of the combustor which entails 
flow and turbulence measurements in nonreacting flows and 
detailed measurements of velocity, temperature and composi- 
tion as well as flow visualization in combusting flows, is 
key to developing new and viable combustor designs for prac- 
tical applications. 

1,3*1 The Swirl Combustor 

The combustor under investigation (Figure 1-1) is a lean, 
premixed combustor in which the flame is stabilized by a 
swirl-generated recirculation zone. The combustor consists 
of two co-axial swirling jets confined in a long cylindrical 
test section of 10,2 cm I.D. The inner jet, ^.8 cm in diame- 
ter, consists of premixed fuel and air. Swirl in the inner 
jet is generated by tangential air injection well upstream of 
the combustor inlet. The outer jet flow is air only and has 
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an outer diameter (D^) of 10,2 cm. Swirl in the outer jet, 
which is generated by varying-angle vanes upstream of the 
combustor inlet, can be either co-swirl (outer and inner 
swirls in the same direction) or counter-swirl (outer and 
inner swirls in the opposite directions). The fuel, commer- 
cial grade methane or propane, is injected radially from the 
inner Jet centerline five combustor diameters (5D^) upstream 
of the inlet. Ignition is accomplished by a long reach spark 
plug, A more detailed account of the combustor can be found 
in Chapter 2. 

The swirl number as defined in Section 1,1 is non-dimen- 
sional ized by the outer radius of the jet. The swirl number 
for the inner or the outer flow is then given by 


S 


*^uvrdr 


U/" .. 2 


2 n r 7 u*" r dr 


(1.6) 


where r^ and r^ are the inner and outer radii of the ‘jet 
under consideration, u and v are the time-mean axial and tan- 
gential velocities respectively. The swirl number is posi- 
tive as defined, A negative sign is used for the outer swirl 
number under counter-swirl conditions. 


The principal operating variables of the combustor are the 
inner and outer swirl numbers, the inner and outer jet flow 
rates and the fuel flow rate. The axial velocity ratio is 
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the ratio of the volumetrie mean axial velocity of the inner 
jet to that of the outer jet# The combustor is operated at 
atmospheric pressure without preheat. 

• 3* 2 Results from Previous Studies 

A series of studies as part of the overall program have 
been made examining different aspects of combustor operation. 
Martin[383, who also built the combustor, and Martin, et 
al.C39] measured stability and blow-off limits for methane in 
the combustor. Their[38,39] results show that the lean 
blow-off limits are very close to the lean inflammability 
limits and are independent of outer swirl. The rich blow-off 
limits, however, are lower than the rich inflammability lim- 
its and show significant variations with outer swirl, the 
limits being higher for higher outer swirl numbers and higher 
axial velocity ratios, thus indicating the importance of 
swirl-generated recirculation zone size and turbulent mixing 
for the rich limits. AlwC, flame instability was encountered 
in rich mixtures and combustion instability is the likely 
cause for blow-off at the rich limits. Their[38,39] visual 
observations are the following: co-swirl cases have long 

cylindrical flames while counter-swirl cases have much 
shorter bubble^like flames. The flames for lower axial vel- 
ocity ratios are smaller in size than flames under high axial 
velocity ratios. They conclude that high rich blow-off lim- 
its and best flame stability are achieved under high counter- 
swirl conditions. 
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Oven and Oven, et al .C33» 40, ^1 3 have made NO-NO , tempera- 

A 

ture and species concentration measurements at various axiai 
and radial locations in the combustor for methane firing. 
Yetter and Yetter and Gouldin[42, 433 have made temperature, 
velocity and composition measurements at the exhaust with 
methane as fuel. These investigators[33, 40, 41 , 42, 433 report 
large concentrations of unburned fuel and CO in the exhaust 
which they attribute to quenching in the inter-jet shear 
layer, but NO emissions are low. Oven, et al.[333 report 

A 

large fractions of NOg/NO^^ in the exhaust, especially in 
counter-swirl cases, which they attribute to quenching in the 
shear layer and turbulent mixing. Yetter ’ s[ M2 3 and Yetter 
and Gouldin * s[M 33 results show poor combustion efficiencies 
in general, though efficiency is significantly improved in 
going from high counter-swirl to high co-swirl and as the 
axial velocity ratio approaches one. These results are dis- 
cussed in greater -detail and compared with results from the 
present work in Chapter 3. 

Lee[443 made Laser Doppler velocity measurements under 
cold flow (without combustion) conditions. Results show a 
recirculation zone for the counter-swirl case, but no such 
zone is seen in the co-swirl case in cold flow. However, 
centerline profile measur ements[M43 under combusting condi- 
tions show that a recirculation zone develops for the co- 
swirl case as well, while little alteration of the recircula- 
tion zone due to combustion is seen in the counter-swirl 
case. Visual observations of the ignition sequence, with 



14 


video photography at 20 frames per second, made by 

Halthore[453 show, in combination with Lee’s measure- 
ments[443, that in the counter-swirl case the flame is stabi- 
lized in front of the recirculation zone whereas in the co- 
swirl case the flame is first stabilized in the wake of the 

ignitor electrodes, a recirculation zone is then formed and 

the flame is stabilized in front of this recirculation zone. 
Halthore[453 also made blow-off measurements with propane in 
order to compare with Martin's results for methaneC383. He 
found that the lean blow-off limits for propane are nearly 
the same as those for methane and are also independent of the 
outer swirl. The rich blow-off limits for propane are how- 
ever slightly higher than those for methane. Beyler[46] and 
Beyler and Gouldin[473 made emission spectroscopy measure- 
ments of the flame for methane firing in order to character- 
ize the flame structure and the mechanism of flame stabiliza- 
tion. The results of his investigation establish that the 
stabilization of combustion occurs in front of the recircula- 
tion zone in a region of low velocity which, in the mean, is 
largely free of reci»*oulation zone influence through heat and 
mass transfer. Beyler also found that the structure of the 
reaction zone in the forward regions of the combustor is the 
same for both co- and counter-swirl conditions but differ- 
ences are seen further downstream. 
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^ ^ Present Work 

The present work is one of the steps in the overall 
research program outlined earlier in this section. In the 
first phase of the present study the combustor is operated 
with prevaporized propane. Temperature, velocity and compo- 
sition measurements including NO-NO^^ are made for varying 
swirl conditions along the radius of the combustor at the 
exhaust plane, six combustor diameters (6D^) doi^mstream of 
the combustor inlet. Temperature measurements are made with 
an uncoated Pt/Pt-10?[Rh thermocouple. A two-holed pressure 
probe is used for velocity measurements. Composition mea- 
surements other than those of oxides of nitrogen are made 
with a gas chromatograph. NO and NO are measured with a 

A 

chemiluminescent analyzer used in conjunction with a NOg to 
NO converter. A cylindrical water-cooled probe is used to 
draw gas samples for composition measurements. From these 
data chemical and thermal efficiencies are calculated for 
various operating conditions. 

Propane is chosen as the fuel for two reasons. First, 
propane is more representative of the commercially used fuels 
than methane, the previously studied fuel. Second, the 
results for propane can be compared with the previous results 
for methane (see Section 1.3.2) to investigate the effect of 
possible differences in the chemical processes on the combus- 
tion of the two fuels in the combustor thereby getting a bet- 
ter insight into the operating characteristics of the 
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combustor. In the second phase of the study cold flow 
measurements (without combustion) are made for propane in 
order to determine the distrubution of the fuel in the 
exhaust plane as well as at the inlet to the combustor. Cold 
flow measurements are made for methane fuel also. In addi- 
tion, the measurements made in the first phase of the present 
study are repeated using methane fuel. These results are 
compared with previous results for methane (Section 1.3*2) as 
well as the results for propane obtained in the first phase 
of the present study. 

Chapter 2 of this thesis contains the details of the 
experimental set-up and the measurement techniques. Results 
are presented and discussed in Chapter 3* Chapter 3 also 
contains details of the calculation of combustor efficien- 
cies. The conclusions from the present study are noted, and 
suggestions for future studies are made in Chapter M. 


CHAPTER 2 


EXPERIMENTAL 

2. 1 Test Facility 

A schematic layout of the test facility is shown in Fig- 
ure 2-1. Ambient air is supplied to the combustor by a sin- 
gle blower, a belt driven Buffalo Forge centrifugal fan, 
model 7E, run by a 25 HP General Electric motor. Air sup- 
plied by the blower is divided into two streams the inner and 
outer flows for the combustor. Downstream of the blower each 
stream passes through a metering station, a settling chamber 
and a variable swirl generator. Butterfly valves in each 
line permits adjustment of the inner and outer flow rates 
which are measured with annubar elements (Ellison Instrument 
Division, model 710, Dietrich Standard Co.). Temperature 
measurements by copper-constantan thermocouples and pressure 
measurements from wall pressure taps are recorded at each 
metering station allowing for density corrections of the mea- 
sured flow rates. Area integrated mass flow rates, obtained 
under nonreacting conditions without swirl, from the combus- 
tor exit velocity and temperature profiles were used to det- 
ermine the calibration constant for each annubar. The accu- 
racy of the calibration is estimated to be - 2 % of the 
instantaneous flow rate. These calibrations were originally 
made by Martin[38] and were repeated by subsequent research- 
ers who have worked on the test rig. The calibrations were 

not repeated in the course of the present work. 
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straightening sections are placed upstream of the metering 
station to eliminate radial and tangential velocity compo- 
nents and to obtain fully developed flow profiles, 

Both inner and outer flows pass through separate settling 
chambers and variable swirl generators (Figure 2-2). Air 
leaves the variable swirl generators and enters the combustor 
through two confined co-axial jets . The inner jet, ^.76 cm 
in diameter, is premixed with fuel and the outer annular jet, 
10.16 cm in diameter, is air only. The combustor test sec- 
tion is 58.3 cm long, 10.16 cm in I.D., and made of quartz 
tubing. The inner diameter of the test section is the 
’combustor diameter' (D^) . Variable swirl in the inner flow 
is obtained by injecting air both tangentially and axially 
into a mixing section approximately five combustor diameters 
upstream of the combustor inlet. Swirl in the outer jet is 
regulated by variable angle vanes with both co- and counter- 
swirl conditions possible. A long reach spark-plug is 
inserted near the centerline close to the inlet for the pur- 
pose of ignition. The exhaust gases exit through a 10.16 cm 
I.D. stainless steel tube section and out through a chimney. 
Measurements of gas composition, velocity and temperature in 
the exhaust were made at a station six combustor diameters 
downstream of the combustor inlet, by means of probes 
inserted radially and supported by the stainless steel 
exhaust tube. 


The fuel, either commercial grade propane or methane, is 
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injected radially from the inner jet centerline 5D^ upstream 
of the combustor inlet through eight orifices symmetrically 
positioned around the circumference of a 0.1|5 cm O.D. stain- 
less tube. Fuel metering is accomplished by a Brooks (Type 
1 1 1 0-08H2G 1 A ) rotameter calibrated for propane at 
4.8x10^ N/m^ (69*7 psia) and 21^C. For other metering condi- 
tions density corrections are made using the temperature 
obtained with a Cmega Engineering (Model 199) iron-constantan 
thermocouple and pressure obtained with a bourdon-tube (Maxi- 
safe A.I.S.I 316) pressure gauge. Two propane cylinders par- 
tially filled with liquid propane supply the fuel. The two 
cylinders are placed in a water bathj whose temperature can 
be controlled by a flow of suitable amounts of hot and cold 
water from the building water lines. The water bath supplies 
the heat for vaporization of the fuel. (See Reference 42 for 
methane metering.) 

The combustor, which is operated at atmospheric pressure 
without preheat, has a maximum overall inlet velocity of 
45 m/s. A more detailed account of the combustor design, 
construction and instrumentation has been given by 
Martin[38]. Equations for determining the combustor operat- 
ing conditions are listed in Appendix A. 

2. 2 Composition Measurements 

Exhaust gases are sampled from radial locations in the 

exhaust plane (6D^ downstream of the inlet) and analyzed for 

nitric oxide (NO), total oxides of nitrogen (NO v , 

^ x) , carbon 
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dioxide (COg), carbon monoxide (CO), unburnt propane (C^Hg), 
lower hydrocarbons (e.g., ethane (CgHg), ethylene (CgHij) and 
methane (CH^)), oxygen (Og) and nitrogen (Ng). 

2.2.1 Sampling Probe 

A water-cooled cylindrical probe spanning the diameter of 
the combustor is used to sample the exhaust gases in this 
experiment. A transverse cylindrical probe is felt to have 
certain advantages over the conventional L-shaped probe. For 
a cylindrical probe, in flows with very low radial veloci- 
ties, sample hole alignment to the flow is accomplished by 
simply rotating the probe to balance pressures at two taps 
located at equal angular displacements on either side of the 
sampling hole on the surface of the cylinder. Also the dis- 
trubance to the flow is more symmetric with a cylinder span- 
ning the diameter than an L-shaped probe inserted into the 
flow. Finally, the cooling passage design is simplified, 
resulting in a smaller overall probe diameter since cooling 
water flows in one direction only. 

OvenC^O] considered various materials and designs for the 
probe and found that the probe shown in Figure 2-3 was best 
suited for sampling in the exhaust plane. Considerations in 
the choice were coolness of the probe to prevent catalytic 
reactions in the hot center piece (i.e. the thermal conduc- 
tivity of the center piece material) and minimum probe inter- 
ference to the flow. 



21 


The same probe shown in Figure 2-3 is used in the present 
work. Two pressure taps, a sampling hole and a cooling water 
hole are drilled into a solid cylindrical brass piece of 
7.9 mm diameter. Looking at the cylindrical cross-section, 
one sees that the pressure taps are located 59° on either 
side of the sampling hole. Pieces of stainless steel tubing, 
silver-brazed to this center piece, form the pressure and 
sample lines. Stainless steel tubing of 7*9 mm outer diame- 
ter forms the outer jacket of the probe* providing a cooling 
water flow passage as well as structural support. The sample 
line of the probe is connected to the sample train which 
draws the sample for analysis. The pressure lines are con- 
nected to either ends of a micromanometer enabling pressure 
measurements for velocity calculations as well for alignment 
of the probe to the flow direction. 

The need for rapidly quenching chemical reactions in the 
gas sample when once inside the probe, to preserve sample 
composition, is widely r ecognized [ 30 , ^2, ^8 , ^9 ] . Hence the 
probe is cooled by a water flow through it. Cooling water is 
supplied to the probe at room temperature. 

GouldinC50] has shown that low velocity sampling or isoki- 
netic sampling are not satisfactory techniques for sampling 
in regions of large density and velocity fluctuations such 
as, for example, turbulent reaction zones, since the effect 
of the above turbulent fluctuations is to cause the probe to 
draw a sample which is not truly representative of the mean 
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composition in the flow. He suggests either choking the flow 
at the sampling hole or sampling with a large pressure drop 
across the sampling orifice* In the present work a bellows 
pump provides a large pressure drop for sampling bub the flow 
is not choked in general. Moreover, the sampling is done in 
the exhaust plane wliere density fluctuations are much lower 
than in the reaction zone. 

Alignment tests were performed by Oven[ilO] to determine 
the effect of probe (sample hole) misalignment with the flow 
direction on the composition measurements. His results 
showed that probe misalignment upto 60 degrees had little 
effect on the composition measurements in the present combus- 
tor. In the present experiment, however, the probe is 
aligned with the direction of the flow (see Section 2.3 for 
the method used to determine the flow direction). 

The probe is supported horizantally across the diameter of 
the combustor by two brass sleeves, fitted diametrically 
opposite each other, to the walls of the combustor. A scale 
attached on the outside, to the combustor wall, and a dial 
protractor attached to the cylindrical probe serve to posi- 
tion the sample hole and pressure taps at desired radial 
locations and to determine the angle of rotation (about its 
own axis) of the probe. 

2.2.2 Sample Train 


Figure 2-^ shows the sample transfer system. As with the 
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probe it is necessary to prevent any reactions in the sample 
transfer lines in order to preserve sample composition. The 
same material and temperature considerations pertinent to the 
probe apply to the sample lines. All lines are 6.4 mm O.D. 
stainless or polyethylene tubes with stainless steel and 
brass fittings. 

• 

The main sample line is connected to the probe sample line 
and thus draws the sample into the transfer system. The main 
line is then divided into branches — (a) the NO analyzer 

A 

line and (b) the syringe sampling line. In line (a) a stain- 
less steel diaphragm pump (Metal Bellows MB-41) draws off 
0.15 1/min. of sample and pressurizes it to atmospheric pres- 
sure to pass through the NO converter. From the converter 

A 

the sample is bled across a needle valve to the chemilumines- 
cence reaction chamber which is maintained at 5 torr by a 
vacuum pump. The main line and the sample lines upto the 
reaction chamber are resistance heated to 50°C by nichrome 
wire, to prevent removal of NO 2 by water condensation in the 
lines. It is expected that potential reactions in the sample 
lines are avoided since stainless steel tubes are nonreactive 
with respect to oxides of nitrogen at temperatures less than 
lOO^C [51]. The lines are electrically insulated with asbes- 
tos material, wrapped with nichrome wire, and wrapped with an 
additional layer of asbestos insulation. 

Sample gas is drawn into the syringe sampling line (b) by 
another metal bellows pump (MB-41). Since water gives long 



bailing peaks on most chromatographle columns, the sample 
line (b) posses through two ice baths to condense out the 
water. Initially a mixture of dry ice and acetone was used 
in the second bath, but it was found that an ordinary ice 
bath was able to condense out all the water in this experi- 
ment, The syringe sampling port is a septum fitted kajon 
fitting in the sampling line. 

2.2,3 NO, NOj^ Measurements 

NO and NO are measured with a chemiluminescent analyzer 

A 

constructed at Cornell University in conjunction with a 
Thermo Electron (N-CV-1 526-36) NO^^ converter (Figure 2-^). 
The NO converter, a stainless steel tube heated to 1023°K, 

A 

converts NO^ to NO. The chemiluminescent analyzer is sensi- 
tive to NO only. In combination with the converter the total 
NO can be measured; when the converter is bypassed, NO alone 

A 

is detected. The difference between the above two measure- 
ments gives the amount of NOg in the sample. A brief but 
comprehensive description of the principle and operation of 
the chemiluminiscent analyzer can be found in Reference 52. 

In the present study the reaction chamber is maintained at 
5 torr by means of a vacuum pump, with approximately 
0.15 1/min sample flow and 0.03 1/min ozone/oxygen flow. The 
ozone flow is obtained by passing oxygen through a 15 KV 
alternating electric field in an ozonator constructed at Cor- 
nell. Though a measurement of the amount of ozone produced 
was not made in this study, it is expected that a mixture of 
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about 2.0 percent ozone in oxygen leaves the ozonator accord- 
ing to tests conducted by HomanC033. The oxygen pressure in 
the ozonator is about 1,13x10^ N/m^ (16.^ psia) . The output 
current from an uncooled EMI photomultiplier tube operated at 
900 volts is converted into a voltage drop across a 
5,6 Megohm resistance and read out on a Dana digital voltme- 
ter. The NOj^ converter is operated at atmospheric pressure. 
It was observed that the output signal from the photomulti- 
plier tube is very sensitive to the pressures in the ozonator 
and reaction chamber and the sample and oxygen flow rates. 
Hence care is taken to keep these readings steady and cons- 
tant throughout a test run. 

As mentioned earlier, all lines in the chemiluminescence 
system and the sample transfer lines are maintained at 50°C 
(323°K) by heating to prevent removal of NOg by water conden- 
sation in the lines. Significant reduction in measured NOg 
has been observed when water vapor is allowed to condense in 
the lines[03, 5^3. Tuttle, et al.[553 have reported that par- 
ticulate matter, especially soot, trapped in sample filters 
affect NOg measurements since NOg is absorbed or reduced by 
the carbon and hence larger analysis times are required 
before accurate NOg readings are obtained. In the present 
case, however, no particulates are evident in the exhaust 
because of premixing and the nature of fuel used. Even so, 
about two minutes are provided for each measurement to allow 
the detector to reach its full response and to account for 
any effects of particulates which could have been trapped in 
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the system during use of the onalyzi^r in other experiments. 

It is now well known that the presence of CO, COg, hydro- 
carbons and Hg cause interference in NOj^ measurements by 
reduction reactions in which NOg and NO are reduced to Ng and 
perhaps 56, 573. Interference depends on converter 

history and conditioning, oxygen level in the sample or ove- 
rall equivalence ratio and converter temperature. Stu- 
diesC55, 56, 573 indicate that Interference is serious only 
when the levels of the above mentioned interfering compounds 
are high, the fuel-air mixture is rich or when the converter 
is not conditioned. The results of Sigsby, et al.[563 and 
Breithenbach and Shelef[583 which are also listed in Refer- 
ence 55 suggest that as long as the sum of CO and HC concen- 
trations in the sample is less than the oxygen level and the 
converter is conditioned at close to operating temperature by 
flowing air through it for about two hours, accurate N0„ mea- 
surements can be obtained with the converter operating at 
1023°K. Sigsby, et al.[563 suggest that preconditioning will 
oxidize the metallic surface of the converter and prevent the 
exposure of any reducing agents on the surface to the sample 
stream. In the present experiment Hg is not present in the 
exhaust gas, combined levels of CO and HC are far less than 
oxygen levels in the sample and the sample gas is oxidizing 
(oxygen rich) in character. Hence all the above conditions 
for reliable and accurate NO measurements are satisfied. The 

A 

converter is conditioned at 1023°K with air flow for at least 
one hour prior to every test run and operated at 1023°K, In 
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addition the entire system including the converter is purged 
with air between every data point. In this way oxidation of 
the metal surface of the converter is ensured. 

Some attenuation in the detector output is caused by col- 
lisional deactivation of NOg by Og, N 2 > COg and HgO. Deacti- 
vation by Ng is accounted for in the detector calibration. 
To correct for the remaining species, a quenching efficiency 
is calculated for each species based on its respective 
cross-sGction for collisional deactivation of NOg and its 
concentration in the sample. Based on the method described 
in Reference 59» calculations by OvenC^O] show that the com- 
bined effects of Og, COg, and HgO attenuate the detector out- 
put by less than 3 %- It is expected that the results of the 
calculations will be no different for the present experiment 
due to the similarity in the sample gas composition with 
Oven’s experiment. Moreover, the sample is predominantly Ng 
whose effects of collisional deactivation are accounted for 

in the calibration. Since the effects of collisional deacti- 

^ * 
vation of MOg due to other species is expected to be minimal 

based on Oven’s calculations, no corrections are made to the 

NO-NO data presented in this thesis. 

A 

The analyzer is calibrated before and after every run. 
The linearity of the instrument is well accepted. Hence a 
straight line calibration curve is drawn between the voltage 
corresponding to the photomultiplier dark current and that 
corresponding to ^23 ppm of NO in Ng calibration gas. The 
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analyzer can Indicobe the ooncenfcrnbion bo wibhin 0.5 ppm 
based on bhe noise in bhe volbage signal. 

2 . 2 . ^1 Analysis by Gas Chromabography 

The objecbive of bhe analysis is bo obbain bhe relabive 
amounbs, in a dry mixbure of oonibusbion producbs of propane, 
of oxygen, nibrogen, carbon monoxide, carbon dioxide, unburnb 
propane and obher possible hydrocarbons like ebhane ,ebhylene 
and mebhane, in a simple way preferably wibh a single injec- 
bion of sample and a shorb analysis time. Samples are brans- 
ferred from bhe sample line bo bhe chromatograph by syringe. 

Sample Transfer Via Syringes 

The sampling port in the sample transfer line is equipped 
with a septum allowing insertion of a syringe needle for sam- 
ple extraction. Hamilton 1005 5 ml gas bight syringes are 

used to remove and store samples before injecting them into 
the gas chromatograph for analysis. Syringes are flushed 
while inside bhe sampling line by repeatedly drawing • and 
ejecting bhe gas sample, before finally drawing bhe sample 
for analysis. The septum in bhe sampling port maintains 
air-tightness while flushing bhe syringe and collecting the 
sample. About 2.5 ml of gas sample is slov^ly drawn into bhe 
syringe over a period of about 10 seconds. The pressure in 
the sample line at the sampling port is just above atmos- 
pheric. Hence when the syringe is taken out of bhe sampling 
port, outside air cannot enter the syringe due bo the 
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slightly higher pressure inside the syringe. The syringe is 
then immediately inserted into a specially made stand such 
that the needle tip pierces a septum, at the bottom of the 
stand which blocks the needle hole and prevents air diffusion 
into the syringe. 

A maximum of 9 samples are taken during any one test run. 
All syringe sample’s are analyzed within 2.5 hours of sampling 
to minimize errors due to leakage. Samples stored for longer 
periods of time showed, upon analysis, significant air leak- 
age into the syringe. Oppegaard [ 60 ] has identified the vari- 
ous problems facing a syringe injection with septum and sug- 
gested ways to overcome them. The procedures followed in the 
present analysis concur with most of his suggestions, except 
that the new injector design suggested by him is not used. 

Selection of Column and Technique 

The problems facing a simplified analysis involving a sin- 
gle injection are the following. Firstly, there is no single 
column material which allows complete separation of all the 
aforementioned gases in the sample, without requiring a long 
time for analysis (50 minutes or more). Secondly, a few 
column materials irreversibly adsorb some gases. For exam- 
ple, molecular sieve, which is the most efficient column to 
obtain nitrogen/oxygen split, irreversibly adsorbs carbon 
dioxide,. Porous polymers and silica gel which allow separa- 
tion of carbon dioxide are incapable of separating oxygen 
from nitrogen and of separating carbon monoxide from the com- 
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posite air peak at room temperature or above. 

To overcome the problem of separating carbon dioxide and 
carbon monoxide as well as oxygen and nitrogen, Alltech Asso~ 
dates provides a concentric column (CTR column) [61] which is 
a column within a column and thus permits the use of two 
different packings in parallel. Carbon dioxide is separated 
by the inner poropak column while oxygen , nitrogen and carbon 
monoxide are separated (in that order) by the outer molecular 
sieve column. The separation of any mathane in the sample is 
also accomplished by the molecular sieve column with the 
methane peak occuring between the nitrogen and carbon monox- 
ide peaks. The entire analysis is isothermal, takes about 6 
minutes and is accomplished with a single sample injection. 
Tiie CTR column, however, is unsuitable for the samples of 
interest in the present analysis. Firstly, ethylene and 
ethane eluted from the inner poropak column overlap the oxy- 
gen and nitrogen eluted from the molecular sieve column. 
Secondly, though propane can be eluted from the poropak 
within 10 minutes by temperature programming, it takes nearly 
80 minutes to elute from molecular sieve. Thus a subsequent 
injection of sample is delayed for 80 minutes. 

The possibility of two separate injections on two diffe- 
rent columns, one a molecular sieve and the other a poropak, 
was also investigated. In this scheme the molecular sieve 
would separate 0^, Ng, CO and CH^, while the second injection 
into poropak enables the separation of CH||, CO^, OgHg 
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and CgHg. This scheme also poses problems. The propane 
injected into the molecular sieve column is retained for a 
long time during which a second injection is not possible. 
In addition, to obtain concentrations on a relative percent 
basis with schemes involving two separate injections, either 
exactly equal amounts of sample should be injected with each 
of the injections or there should be at least one common peak 
in the two chromatograms whose relative areas can then be 
used to relate the areas of the other peaks in the chromato- 
grams. In the above scheme, the only common peak is that of 
methane which may or may not be present in the actual samples 
obtained from the combustor. Even if present, it may be in 
very small amounts, whereas it is desirable that the common 
peak be a major component in the sample. Also, injecting 
precisely the same amount in both the injections using syr- 
inges is difficult. 

Kaiser[62, 63 ] synthesized a new material called the car- 
bon molecular sieve, for use in gas chromatography. The use 
of carbon molecular sieve for analyses of various gas samples 
have been demonstrated by Zlatkis, et al.[64], Bollman and 
Mortimore[65 ] and others. Unfortunately, none of the ana- 
lyses produced the oxygen/nitrogen split and at the same time 
gave elution of propane in a short time. Bollman and korti- 
more[65] used carbon molecular sieve marketed by Supelco Inc. 

( Bellefonte , PA ) under the name of carbosieve B[66]. Supelco 
introduced a new type of carbosieve called carbo- 
sieve S[67,68] in 1976. Carbosieve S is spherical in shape 
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and has increased surface area. It is best suited for the 
separation of permanent and sulphur gases. Carbosieve S 
gives separations similar to carbosieve B and hence Carbo- 
sieve S alone cannot solve all the problems posed in the pre- 
sent analysis. Using a carbosieve S column and temperature 
programming upto 230°C (503°K) it is possible to obtain com- 
plete separation of all the gases in the sample of interest. 

But in order to get a good oxygen/nitrogen split the column 
must be long and the initial temperature must be close to 
room temperature. These conditions increase the time of elu- 
tion of propane to approximately ^10 minutes resulting in a 
total cycle time (including cooling of the oven) of ^5 to 50 
minutes . 

It is crucial to keep the analysis time per sample, short. 

The combustor is allowed to run for a while to reach a steady 

state of operation. The samples from about 9 radial loca- 
tions are then collected and stored in syringes, over a per- 
iod of i|5 minutes to 1 hour. It is desirable not to run the 
combustor for longer periods of time in order to conserve 
fuel. If each sample analysis were to take 50 minutes the 
total time for analysis of all the samples collected in a 
single run would be 7 to 9 hours. Substantial leakage of 
ambient air into the syringes during this period of time 
would cause unacceptale contamination of the last few sam- 
ples. Though it is possible to store samples in evacuated 
bottles or sampling bulbs, which makes the sampling system 
more complicated, it is nevertheless desirable to have a 
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short anoXysis time. 

To overcome all these difficulties a scheme using two 
columns in series, with a ^l-port switching valve is used. 
After a careful study of the properties of various columns, 
it was found that columns of poropak Q and carbosieve S would 
be most suitable for the purpose. 

Figure 2-5 shows a schematic of the column arrangement and 
the flow path. The ^-port valve switches between the two 
flow arrangements, ’A' and *B'. In position ’A' the two 

columns are in series and the elutants from the poropak Q go 
through the carbosieve S column. In position the carbo- 

sieve column is bypassed and the elutants from poropak go 
directly to the detector. 


The gas chromatograph conditions are as follows. 


Column 1 
Column 2 

Oven temperature 

Carrier gas 
Flow rate 
Detector 

Sample size 


Poropak Q, 1/8"x6 ft., SS 

Carbosieve S, 1/8*'x8 ft., 100/120 mesh, SS 

50®C for ^ minutes and programmed 

at 25°C/min. to 175*^C 

Helium 

^7 ml/min. 

Thermal Conductivity, 270‘^C, 160 MA, 

0.5 mV range 
0.5 ml . 


Figure 2-6 shows a chromatogram obtained under the above 
conditions for a sample withdrawn from the combustor. The 
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time axis between the peaks in the chromatogram is compressed 
to accomodate all the peaks. 

The sample is injected with the 4-port valve in position 
'A*. A composite peak of air, carbon monoxide, methane (if 
present) and carbon dioxide elutes within 0.8 min. from the 
poropak and passes on into the carbosieve column. The valve 
is switched to position ’B' at the end of 2 minutes. Ethy- 
lene and ethane, which elute in that order at 2.5 minutes 
from poropak, pass directly to the detector. (Ethylene and 
ethane are not present in the sample and hence are not seen 
in Figure 2-6). 

The valve is switched back to position ’A* immediately 
after the elution of ethane (i.e. after approximately 3 
minutes) . Oxygen elutes from the carbosieve column at about 
3.5 minutes (the exact time depends on how long the valve is 
kept in position 'B’), followed by nitrogen and carbon monox- 
ide. The valve is again switched at the end of 6 minutes to 
position ’B'. Propane eluting from the poropak passes 
directly to the detector and can be detected at about 6.8 
minutes. The valve is then switched back to position ’A'. 
Methane elutes from the carbosieve at 9.5 minutes followed by 
carbon dioxide at 11.2 minutes. The column oven temperature 
is then reset to the initial temperature, 50°C. In this way 
the entire sample analysis including propane, ethane and 
ethylene can be done within 12 minutes. In addition, the 
columns are ready for the next sample as soon as they reset 
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to the initial temperature. Hence the total cycle time is 
only 15 to 16 minutes. 


Data Reduction and Error Estimation 


The gas chromatograph used for the composition analysis is 
a Varian Model 3700. The detector output is recorded on a 
Varian Model 9176-02 strip chart recorder which has a 
built-in triangular trace electronic integrator. The number 
of transverse line crossings of the integrator trace is pro- 
portional to the area of the peak. The operating conditions 
and the columns used are described in the previous subsection 
on technique. Area data of the peaks are reduced by a rela- 
tive percent method to obtain concentrations as mole frac- 
tions. In order to do this, nitrogen, being the major spe- 
cies in the sample, is chosen as reference and a response 
factor is obtained for each of the other species by calibra- 
tion using calibration gases. Calibration gases used are 
mixtures of the gas species of interest, whose concentrations 
are known to within *255 accuracy. The calibration gases were 
obtained from Alltech Associates and their compositions are 
listed in Table 2.1. Response factors (b^) are defined as 


ResponseC peak area) per unit concentration of 

bi= 

Response per unit concentration of species i 

which gives. 



^N2 ^N2 ^ ^cN2 
A. a^ / X^. 


( 2 . 1 ) 


( 2 . 2 ) 


where : 
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A = the area of the peak 

a = detector attenuationC set on the gas chromatograph) 
Xq = mole fraction of the species in the calibration 
gas sample. 


It can be seen from the above expression that response 
factors can be determined for each species in a calibration 
chromatogram indepedent of the the other species. The res- 
ponse factor for nitrogen is 1 by definition. With these 
response factors the mole fractions(xp of the various spe- 
cies in a sample gas can be determined by, 


= 


^ 

f "'j "J *’J 


(2.3) 


In this method the area of each species is converted into an 
quivalent area of a nitrogen peak of the same concentration 
as that of the species. The numerator in Equation 2.3 repre- 
sents an equivalent nitrogen peak area for species i. The 
denominator is the sum of all such areas for each species 
including nitrogen. Hence the resulting concentration will 
be in mole fraction units. 


The advantage of the relative percent method is that it is 
independent of the sample volume injected. A disadvantage is 
that an error in measuring the area of one component peak 
will propogate through out the numbers for the other compo- 
nents as v/ell. Also, since the total of the concentrations 
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must always add up to 100 percent, unmeasured species con 
introduce error if their concentrations are significant. 
Results presented in this study should be free of errors of 
this type since all the species present in any significant 
amount are measured. 

A number of calibration injections are done for each of 
the calibration gas mixtures and response factors are calcu- 
lated for every calibration injection. These data were sta- 
tistically analyzed to obtain standard deviations for the 
response factors. The mean of the response factors for each 
of the species is used in calculations using Equation 2,3- 
Error analysis of Equation 2.3 shows that contributions to 
the error in come mainly from the error in the response 

factor for the species i and the uncertainity , arising from 
the detector noise and the least count of the integrator 
markings on the chart paper, in determining the area of the 
peak for the species i. Neglecting the higher order contri- 
butions from cross terms from other species in the sample, 
the fractional error in is given by 

where and are fractional errors in bj^ and 

Aj. respectively. It is seen from Equation 2.^ that the 
actual error in X^ depends on the actual area of the species. 
Representative concentrations of each of the species were 
chosen to calculate typical values for errors in Xj. . Table 
2.2 lists the values of the response factors and estimated 
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errors in calculated species concentrations, Obviously, the 
errors will be slightly smaller for species concentrations 
larger than those listed in the table and slightly larger for 
lower species concentrations. It should be noted that the 
estimated errors do not take into account the possible air 
leakage into the syringes and the sampling lines which could 
be s serious problem if proper care is not taken to prevent 
leaks. Response factors for ethane and ethylene are not 
listed in Table 2.2 since those species are not encountered 
in the exhaust gases. 

Based on the noise levels, the lowest attenuation and the 
detector response bo the various species, the minimum detec- 
table amounts were calculated to be 40 ppm for CO 2 and aprox- 
imately 70 ppm for the other species. 

Injection technique and consistency in injection are found 
to be critical for good results. Very fast or very slow 
injection results in improper peak shapes and erronous data. 

2. 3 Velocity Measurements 

The pressure taps in the sampling probe are used for vel- 
ocity measurements also. The taps, 0.8 mm in diameter, are 
located 118 degrees apart in the same axial plane. They are 
located symmetrically about the sampling hole which is at an 
axial distance of 8 mm from the plane of the pressure taps. 
The pressure lines from the probe are connected to the two 
sides of a micromanometer. 
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The reference for determining the flow direction (0) is 
token to be the plane containing the flow direction on the 
centerline so that OsO on the centerline. The axisymmetry of 
the 0 values measured from this reference (see Section 3*2) 
supports the choice of this reference. The direction (6) of 
the mean total velocity (neglecting radial components) with 
respect to the reference plane is located by rotating the 
probe about its axis, with the two pressure taps facing 
upstream towards the flow. V/hen the pressures in the two 
taps are equal, as indicated by a null balance on the micro- 
manometer, the local flow direction is aligned with the 
bisector of the angle between the pressure taps. The angle 
made by the bisector with the reference plane gives , which 

is read from the circular dial protractor attached to the 
probe (see Figure 2-7). One of the pressure lines is then 
disconnected from the micromanometer and the pressure at the 
other pressure tap, which corresponds to the pressure on 

the surface of the probe at a 59° angle from the flow direc- 
tion, is read. The probe is then rotated by 59° (see Fig- 
ure 2-7) to align the pressure tap to the local mean flow and 
the total pressure is measured. The difference between 
the two readings gives CP-^p-P^g). Readings are taken along 
the radius of the combustor at 5 mm inter 'als. 

Velocity measurements in turbulent swirling flows using a 
transverse cylindrial probe is affected by stream turbulence 
and swirl in the following manner. Turbulence in the flow 
affects the pressure distribution around the surface of the 





probe in fcwo ways. Firstly, increased turbulence in the flow 
accelerates the transition from a laminar to a turbulent 
boundary layer[C9] and secondly, the separation point for the 
turbulent boundary layer is further downstream on the probe 
surface than for the laminar boundary layer[693* Hence the 
calibration of the probe done in a laminar flow may not be 
valid in a turbulent flow. 

Bilgert^^83 points out the possibility of secondary flow 
along the probe towards the vortex core due to the positive 
mean pressure gradient in the radial direction in swirling 
flows. The effect of the secondary flow is such that the 
velocity measured by the probe at a given radial location 
actually corresponds, in the undisturbed flow, to a farther 
radial location approximately 0.09 of the probe diameter from 
the measurement pointCyO]. Density gradients in reacting 
flows may affect the velocity measurements in a way similar 
to swirl by means of pressure gradients in the flow. 

The diameter of the probe used in this study is small com- 
pared to the diameter of the combustor. Hence swirl and den- 
sity gradients are not likely to cause appreciable distor- 
tions to the flow when the probe is introduced. 

The probe was first calibrated in laminar flow using a 
calibration jet of 2 inch diameter, to obtain a relationship 
between (PT-P59) and (Pj-Pg) where Pg is the static pressure. 
Whether a correction is required and the magnitude of such a 
correction for the turbulent swirling flow in the combustor 



were determinea subsequently. The calibration jet was non~ 
swirling and unoonfined and hence the static pressure was 
atmospheric. Readings of pressures and were taken for 
various velocities of the jet and the relationship between 
(Pf-P^g) and found to be 

(Pt-Ps)i = 0.76075 (PT-P 59 ) + 0.00il2 (2.5) 

in the velocity range of interest, where 
(P^-Pgg) are in inches of water. The subscript 1 denotes 

that the calibration was done in a low turbulent intensity 
(laminar) jet. 

Using the values of P^-P^^ obtained from the probe read- 
ings and Equation 2.5, the values of ^Pt"^S^ 1’ which are the 
dynamic pressures, are calculated. The mean total velocity 
at each radial location is calculated from the dynamic pres- 
sure as follows 

Vj^ = /(2?„BhRgT/Pg) (2.6) 

where 

9 = density of the manometer fluid (water )=1 gm/cc 

g = acceleration due to gravity = 9*81 m/s 
h = dynamic head = (P.p-Pg ) ^xO . 025^ meters 
Rg = Gas Constant for ai'" = 215^.69 (mm-Hg)cm^ /gm°K 
T = local gas temperature in °K 

Pg =local pressure, assumed atmospheric .( see Section 3.3.2). 
Substitution of the above values gives 

V, = ir(1073.9(PT-Ps)iT/P,t'j 


(2.7) 
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where inches of water and is in mm of Hg. 

Since the turbulent intensities in the combustor are high 
(about of the mean flow velocity) there is a need to 

check the validity of the probe calibration (Equation 2,5) in 
the combustor. Laser Doppler Velocity (L.D.V.), described in 
Reference 71, was used to measure velocities in the combustor 
(without the probe in place) and compared with the velocities 
measured by the probe using the calibration in Equation 2.5. 
It was found that the velocities measured by L.D.V. were 
about 207- lower than those obtained by the laminar calibra- 
tion of the probe. Hence the velocity obtained by Equations 
2.5 and 2.7 has to be multiplied by a correction factor, f , 
in order to get the correct velocity, V. It was found that 
f^ was a function of alone. The f^ value- obtained with 
varied outer swirl and at various radial locations under both 
hot and cold flow coiuiitions are plotted gainst in Fig- 
ure 2-8. It is found that f^ changes very slightly with 
and is in the range the range of 0.73 to 0.88 for all operat- 
ing conditions in this study. Since the f values were not 
found to depend on the swirl condition or the temperature 
gradients in the flow, it can be concluded the presence of 
swirl and density gradients in the flow does not affect the 
measurement of velocity using the cylindrical probe. The 
correction factor f^ corrects for any probe blockage effects 
also since L.D.V. measurements are made without the probe in 
place and then compared with the probe measurements to obtain 
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The correcbed mean total velocity (V), and the horizon- 
tal (u) and vertical (v) components of the mean velocity are 
calculated from 

V = V^xf^ , u = V CosO and v = V SinO . (2.8) 

Obviously, calculations of velocities require local temp- 
erature information (Equation 2.7). Measurement of tempera- 
ture is described in the next section. 

Figure 2-9 shows the presi.ure distribution on the surface 
of the cylincrical probe as measured by the two pressure taps 
for a typical flow condition in the combustor. 0 is the 
angle subtended by the bisector between the two pressure taps 
with the flow direction. The Figure shows a symmetric pres- 
sure distribution around the probe. The peaks at +59 and -59 
degrees correspond to the total pressure measured by the two 
pressure taps. Each of the two pressure taps measures the 
same pressure when the bisector is aligned with the flow 
direction (i.e. at 0=0 in the figure) indicating that the 
probe can measure the flow direction accurately. Measure- 
ments of flow direction with L.D.V and with the cylindrical 
probe agree to within half a degree. 

Detailed additional measurements in cold flow comparing 
cylindrical probe velocity measurements with L.D.V. and mea- 
surements with a 5-hole pitot probe previously tested in cold 
swirling flows[72,733 show that while perturbations to the 
flow due to either of the probes are negligible in counter- 
swirl, in co-swirl, under cold flow conditions, both probes 


alter the flow drastically. The velocity profiles measured 
by the two probes and L.D.V (with the probe in place) resem- 
ble those in front of a recirculation zone with strong decel- 
eration on the centerline. The reason for the drastic flow 
alteration is the following. For the co-swirl case in cold 
flow there is no recirculation zone formed near the inlet and 
flow is on the verge of vortex breakdown. Therefore, intro- 
duction of any probe in the flow can cause a drastic flow 
alteration resulting in the formation of a recirculation 
zone . 

?.£ Temperatu re Measurements 

An uncoated fine wire Pt/Pt-10|SRh thermocouple probe is 
used to measure temperatures at radial locations in the 
exhaust plane (Figure 2-10). The thermocouple bead was 
formed by butt welding 0. 12? mm diameter Platinum and Plati- 
num- 10/b Rhodium wires. The thermocouple lead arms are made 
of 0.25^ mm diameter wires for rigidity and support. Mean 
thermocouple e.m.f is obtaind by putting the thermocouple 
output through a TSI averaging circuit. A time constant of 
10 seconds is used and each measurement is averaged for 
approximately 1 minute. Averaging is required due to large 
temperature fluctuations (of the order of 100°K), especially 
in the hot central core of the combustor exhaust. 

The possibility of catalytic reactions on the thermocouple 
junction was investigated by comparing an SiO^ coated thermo- 
couple and an uncoated thermocouple. A difference of about 



5^ was found around 1300°K. Ah higher temperatures and 
higher velocities in the core it was found that the coating 
could not last long enough to obtain any meaningful data. 
However, since the observed difference is small, an uncoatod 
thermocouple is used throughout the experiment. 

The measured temperatures can be corrected for thermocou- 
ple radiation losses to the surroundings, by equating the 
heat transferred from the gas to the probe to that lost by 
the probe due to radiation. The correction is[2^l] 

= eo“d(Tj|-T^)/2k (2.9) 

g c c w 

where 6 is the emissivity of the wire material, o~ the Steph- 
en-Boltzmann Constant (5.6^^x10 erg cm K sec ’ d the 
diameter of the bead and k the thermal conductivity of the 
gas. The temperature subscripts w,c and g refer to the wall, 
thermocouple and gas, respectively. <=<T (t!|“T,^) is the heat 
loss rate by radiation per unit area of the probe. For a 
small diameter spherical bead such that the Reynolds number 
is much less than 1, the heat transferred by the gas per unit 

area of the bead is given by ( 2k/d) (T_-T_) . 

8 ^ 

The emissivity ,£, is a function of temperature. Though 
values of ^ for platinum are not available in literature for 
the entire range of temperature, values listed in References 
7^ and 75, which are collections of data, suggest a value of 
about 0.1 at 1000°K and about 0.2 at 2Q00°K with almost a 
linear variation with temperature. Tlie value of k also 


depends on the temperature. Since air and nitrogen have 
approximately the same thermal conductivity as well as the 
fact that the exhaust gas is mostly nitrogen, the values of 
thermal conductivity for nitrogen listed in Reference 76 for 
various temperatures were used. With a maximum value of 

£=0.2 for the entire calculation and with values of k from 
Reference 76 the temperature corrections were found to be 
89.5®K for T =1923°K and 0.6°K for T =n60°K with T, =3^0°K. 
The corrections are quite small upto about BOO'^K. Hence the 
corrections do alter the temperatures but only in the hot 
central core of the combustor. These corrections did not 

significantly alter the velocities nor the calculated thermal 
efficiency. The velocities on the centerline increases by 
approximately 2?, and the thermal efficiency increases by 
less than 2 %. The reason for this is that the temperature 
correction affects only a small central region and does not 
contribute greatly to altering the efficiency. The tempera- 
ture profiles presented in this work have not been corrected 
for radiation losses unless otherwise noted. 

A method for obtaining corrections for conductions losses 
has been presented by Frock, et al.[77]* The expression is 

(T„-T^) = (T^-T ,)/(Cosh 2a/h/dA ) (2.10) 

g c gw 

where h is the film coefficient of heat transfer for a wire 
of diameter d, A is the thermal conductivity of the wire and 
a is the depth of immersion from the junction to the probe 
supporting tube. The thermocouple used in the present study 
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has a small wire diameter and a long immersion depth of 
1 . 3 cm making the conduction losses negligible. 

2. 5 Experimental Procedure and Test Conditions 

Exhaust composition, velocity and temperature data are 
taken for varying outer swirl, keeping the equivalence ratio, 
inner swirl, axial velocity ratio and overall axial velocity 
at the inlet more or less constant for all the runs. It was 
intended to study variations due to outer swirl since it was 
believed, based on previous results[41 , , that the outer 

swirl was the single most important parameter affecting the 
operation of the combustor. The actual test conditions are 
shown in Table 2.3* 

The apparatus is allowed to warm up and reach steady state 
by running the blower for about half an hour before igniting 
the combustor. After ignition the combustor is allowed to 
run for about 10 minutes to reach steady state operation. 
Measurements are then taken along a radius at intervals of 
5 mm starting from 1 mm away from the wall. The pressure 
taps are first moved to the desired radial location, the flow 
direction determined and the pressure readings taken. The 
sample hole is then moved to the radial location and aligned 
to the flow. After allowing about a minute for the sample 
line to flush completely NO and then NO measurements are 
taken. A sample is collected in a syringe at the sample port 
for gas chromatographic analysis. The pressure taps in the 
probe are then moved to the next radial location and the 




entire procedure is repeated. In the intervals between tak- 
ing data, the NO^ converter is flushed with ambient air 
through a side valve to preserve the oxide coating on the 
converter walls. After traversing the entire radius the 
combustor is turned off to remove the sample probe and intro- 
duce the thermocouple probe. The combustor is re-ignited and 
temperature measurements are made at the same radial loca- 
tions used for velocity measurements. Either probe can be 
positioned at the desired radial location to better than *0.5 
mm accuracy. Gas chromatographic analysis of the sample is 
not made at all the radial locations to keep the total analy- 
sis time short, bub is restricted bo 10 mm intervals and in 
regions of large concentration gradients to 5 mm intervals. 

It was seen after a few runs and analyses of data that 
combustor operation is highly reproducible (to within 3!^ of 
each other as seen by temperature and velocity profiles for 
repeated runs) . Hence it was possible to make individual 
measurements of composition, temperature and velocity from 
separate combustor runs instead of in the same run. This 
reduced the duration for which the samples are stored in the 
syringes . 

The axi-symmetry of temperature, velocity and composition 
in the combustor have been established in previous experi- 
ments by 0ven[40], Yetber[^2] and other authors[44 , 46 ] . 
Hence data in the present study are taken only across half 
the diameter of the combustor and symmetry is assumed in cal- 



culabions. Additional data, taken in the present study, to 
support the assumption of axi-symmetry are presented in Sec- 
tion 3.2, 

In the second phase of experiments, data are taken for 
combustor operation with methane as the fuel. This is done 
for comparison with the propane runs as well as with results 
obtained by yetter[42]. The test conditions for methane are 
also listed in Table 2.3. In addition, composition measure- 
ments are made for cold flow (without combustion) for propane 
at the exhaust plane to study the fuel distribution. Fuel 
distribution at the inlet is also studied for propane and 
methane without combustion. To measure the inlet fuel dis- 
tribution a simple hypodermic stainless steel tube in the 
shape of an L, with the hole facing the flow, is used to 
withdraw samples. 

Swirl numbers are obtained by integration of inlet velco- 
ity profiles measured under nonreacting conditions using a 
five-hole pitot probe. These measurements were originally 
performed by MartinC383 and repeated for calibration checks 
by Yetter[^^2] and Halthore[45] . These measurements were not 
made in the present study. 



CHAPTER 3 


RESULTS AND DISCUSSION 

3 . 1 Introduction 

The results are presented and discussed in two stages; 
results for propane firing (Phase 1) followed by those for 
cold flow and methane firing (Phase 2). The presentation of 
results is proceeded by two other sections. Section 3.2 
deals with the axisyminetry of conditions in the combustor and 
with the reproducibility of results. In Section 3.3 the 
methods for calculating average exhaust composition values 
and efficiencies are given. 

The results for propane firing (Phase 1) are presented in 
Section 3*4. In Section 3.5 the propane results are compared 
with previous results for methane firing[40, 42] . The compar- 
ison shows some disagreements between the two sets of data. 
In an attempt to resolve these disagreements Phase 2 experi- 
ments were performed as follows. Cold flow measurements were 
made and the results of these measurements are presented in 
Section 3.6. However, the results of the cold flow measure- 
ments could not explain the differences. Hence experiments 
were done with methane as the fuel. The methane firing 
results are presented and discussed in Section 3.7. A 
mechanism for the operation of the combustor which can 
explain the previous and the present results is proposed in 
Section 3.8. 
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Composition results are presented on a dry basis in 
Figures (3-1 1 - 3-18,3-25,3-26) and on a wet basis in the 
Tables (3. 2-3. 5 , 3. 7, 3. 8) , 

The folowing terminology will be used in the discussion. 
’Quenching’ refers to the reduction in reaction rates due to 
cooling (thermal quenching) as well as to the freezing or 
stopping of the reactions by the decay and removal of active 
free radicals. Quenching occurs due to the ’dilution’ with 
the cold outer air, of zones where reactions are occuring. 
’Combustion zone’ refers specifically to the region close to 
the combustor inlet (upto about 1.5D^ axial distance from the 
combuster inlet) where most of the combustion is taking 
place. The region downstream of the combustion zone will be 
referred to as the ’post-combustion zone’. The oxidation of 
CO to COg and NO to NOg, which takes several milliseconds and 
is slower by orders of magnitude compared to the fuel 
break-up reactions, occurs in the post-combustion zone. (See 
Section 1.2 in support). 

3.2 Axisymmetry and Reproducibility of Results 

Axisymmetry and reproducibility of results in the combus- 
tor are discussed in Section 2.5. Additional results to 
demonstrate axisymmetry and reproducibility of combustor 
operation are presented in this section. 

Figure 3-1 shows the temperature profile across the entire 
diameter of the combustor. A traverse along each half of the 
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combustor was done from either side of the combustor. The 
temperature values at the same radial distances on either 
side of the combustor centerline differ by less than 3 % of 
each other. Table 3.1 shows the axisymmetry of flow direc- 

tion. (The 6 values on opposite sides of the combustor are 
opposite in sign due to swirl), It can be seen from 
Table 3.1 that the asymmetries in the 0 values are less than 
1 degree. These results justify the assumption of axisymme- 
try at the measurement station. Hence all the results are 
presented for half the diameter of the combustor. 

Measurements in this study were made over a period of 
approximately one year. Hence reproducibility of results is 
important to ensure consistency. Figure 3-2 shows the repro- 
ducibility of a temperature profile. The two profiles were 
measured 4 months apart. The local differences in tempera- 
tures are less than 35S. Table 3.1 also shows reproducibility 
of flow direction measurements to within 1^. 

3. 3 Calculation of Average Values and Efficiencies 

3 . 3.1 Composition on a Wet Basis 

At each radial location the amount of water vapor present 
in the exhaust is determined from a local carbon atom balance 
with the assumption that all H from the burnt fuel is oxi- 
dized to HgO. For every mole of COg or CO in the exhaust, 
4/3 moles of water vapor for propane and 2 moles of water 
vapor for methane are present. Therefore, the amount of 





53 


water vapor in the exhaust for every mole of dry exhaust gas 
is given by 



^’h20 

“ ('^/3)( X'cQg ^’co ^ 

(propane) , 

(3.1a) 


^'h20 

= 2 

( + X’^Q ) 

(methane) , 

(3.1b) 

where 

^’C02 

and 

X’co local mole 

frDotions of CO^ 

and CO 

on a 

dry basis 

. The mole fraction 

of species i on 

a wet 


basis, denoted by > is then given by 


1 + X' 


(3.2) 


H20 


where i refers to 0«, N^, C^H„, CH,, , CO, C0„ or H_0, 

« C 0 0 M 


C. 


3.3*2 Average Values 


Average values of temperature and concentrations in the 
exhaust are obtained on a mass flux weighted basis, Average 
mass fractions are given by 


7 ^ u Y, 2 n r dr 
0 J 1 

YT = , (3.3) 

7 u 2 7^ r dr 


where is the looal time mean mass fraction of species i on 

a wet basis, % is the local density, and (=D^/2) is the 

radius of the combustor. Also where is the 

molecular weight of the species i >nd M is the local molecu- 

lar weight of the exhaust gas mixture. =(PM /RT) where P 

s 



and T ore the local static pressure and temperature 
respectively and R is the universal gas constant, By substi- 
tution Equation 3*3 becomes 


(PM^/RT)uX. (M./Mj2nr dr 
O-' ® iis 

Y. = X. = — . . (3,iO 

M R 

® (PM^/RT)u2nr dr 

0 ^ s 


The static pressure inside the combustor is very nearly 
atmospheric and changes only slightly (by less than 5 cm of 
water column, i.e, 7,4x10^ N/m^) across thr eombustor as seen 
from the pressure readings for velocity measurements. 
Further, the major component of the exhaust gas is air and 
the local molecular weight is approximately equal to that of 
air. For example, at the location where maximum combustion 
products are encountered in this entire study, namely the 
centerline for the high co-swirl propane case (Table 3.2), 
the difference between the molecular weights of the local. gas 
mixture and that of air is only Hence P and are 

assumed constant for the integration in Equation 3.4, which 
then reduces to 



7 (u/T) X. r dr 

Mi 0*^ 


® 7 (u/T) r dr 

0 ’^ 


(3.5) 
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Hence the average (mass flux weighteu) mole fraction in the 
exhaust is calculated by 


(u/T) r dr 



(u/T) r dr 

0 -' 


Based on Equation 3.6 the mass flux weighted average 

exhaust temperature is calculated by 

7 (u/T) T r dr 

0 

T = . (3.7) 

R^ 

7 (u/T) r dr 

Simpson’s rule is used for integration from r=0 to 

r = 5.0 cm, with a step size of 0.5 cm. The velocity at the 

wall (r=R^=5.1 cm) is assumed to be zero. The trapezoidal 

rule for integration is used between r = 5.0 and r = 5.1 cm and 
the value is added to the integral obtained from the Simp- 

son’s rule. An error analysis of the integration shows that 
Simpson’s rule gives the value of the integral to better than 

0.5% accuracy. The reason for this accuracy is that the pro- 

files of the integrands are quite smooth and Simpson’s rule 
which uses a second order polynomial fit is sufficient to 

evaluate the integral accurately. 
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3»3.3 Eff 1 clenoies 


Combustion efficiencies are calculated by three methods. 
In the first method, a chemical efficiency, calculated 
based on the amount of unburnt fuel in the exhaust; 




f 

c 


1 - 



(3.8) 


where is the average mole fraction of fuel (propane or 
methane) in the exhaust calculatd by Equation 3.6 and [F] is 
the average mole fraction of the fuel entering the combustor, 
determined from calibrated rotameter settings. 


The second method is an improvement on the first one. It 
relates chemical efficiency, , to CO in the exhaust as 
well as unburnt fuel since CO not going to CO 2 is lost sensi- 
ble energy. 


n. = 1 - 


^f ^CO^^CO^^f^ 
[F] 


(3.9) 


where q^q is the heat release per mole of CO going to CO^ at 
the inlet temperature , , and is the lower heating value 
per mole of fuel at T^^^. For ease in the actual calculations 
the values of q^Q and q^ at STP ( 67.63 and 484.22 (propane) 
or IS 9.76 (methane) Kcal/gm mole respectively) are used since 
the difference between the inlet temperature and the standard 
temerature is only about 20‘^K and the values of q^^Q and q^ at 
T^n sre not significantly different from those at STP. 
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In the third method the exit sensible enthalpy flux is 
compared with the lower heating value of the fuel to define a 
thermal efficiency, 


'Q. 


/measured increase in sensible enthalpy\ 
\ flux in the working fluid. / 

( lower heating value of the fuel )x m^ 


C3. 10) 


He/tHf^’f f 


where 


R 

H = 9 u (2 Y. h. ) 27Tr dr 

^ qJ ill 


(3.11) 


and 


h. = /c . dT 
1 pi 


T. 
in 


h^ is the increase in sensible enthalpy per unit mass of the 
species i in going from temperature to the local tempera- 
ture T, c . is the specific heat capacity of the species i, 
P 1 

m^ is the fuel mass flow rate into the combustor and q’^ is 
the low heating value per unit mass of the fuel at the inlet 
temperature, The kinetic energy terms u /2 are small 

compared to the species enthalpy changes, h^. , and therefore, 
are neglected. By substituting for and 9 as shown in Sec- 
tion 3 . 3 . 2 , Equation 3.11 simplifies to 


H 


e 


where 

values 



^ (P/RT)u(2 X. / C . dT)2nr dr (3-12) 

0 i 

in 

is the molar heat capacity of the species i. The 
Cp^ as a function of temperature are obtained from 
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sixth order polynomials fibbed by Probhero[78] for various 
species . 

All oalculations were performed with a Digital Equipment 
Corporation PDP-1 1/3^1 computer. A listing of the program 
AVG.FTN is given in Appendix B. The program calculates the 
amount of water, mole fractions on a wet basis in ppm, the 
average values and the chemical efficiencies. The programs 
EFF.FTN and EFFM.FTN which calculate the thermal efficiencies 
for propane and methane firing respectively are also listed 
in Appendix B. 

3 . 3.4 Atom Balance 

As a check on the composition measurements, an atom 
balance was attempted. The amount of water in the exhaust at 
any radial location can be calculated from the dry gas compo- 
sition at that location by two methods. The amounts of water 
calculated by the two methods are compared to check the accu- 
racy of the composition measurements. The first method of 
calculation invokes a carbon atom balance. Based on this 
method the number of moles of water present per mole of dry 
exhaust gas is given by Equations 3.1a and 3.1b. 

The second method is based on an oxygen atom balance and 
the relative percentage of oxygen and nitrogen in air. For 
this method 


X’ 


= 2[(X^2/3.77) - 


H20 


(3. 13) 
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for both propane and methane. Since Equation 3.13 involves a 
difference between two large and nearly equal numbers, small 
errors in composition measurements result in large variations 
in the amounts of water calculated by this method. A statis- 
tical error analysis for Equation 3*13 based on the uncer- 
tainties listed in Table 2.2 yields values of uncertainties, 
for the amount of water calculated, as high as 200^ or more 
especially at large radial distances where small amounts of 
water are present. Similar analyses of Equations 3.1a and 
3.1b yield uncertainties of less than 455. Hence only the 
first method gives reliable values for and is used for 
the calculations. Due to the large uncertainties in the sec- 
ond method a meaningful atom balance cannot be obtvUned. 

Though a check on the composition measurements cannot be 
made at every radial location an overall check on the mea- 
surements can be made by comparing the measured mass fluxes 
of carbon at the exhaust as fuel, CO^ a^^d CO with the inlet 
carbon mass flux. Such a comparison shows that the mass 
fluxes of carbon at the exhaust are lower than at the inlet 
by less than 1055. The lower value for the exit mass flux of 
carbon, calculated from the exhaust gas composition profile, 
may, in part, be due to air leakage into the sample train or 
syringes during composition measurements. 
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Phase 1 : Propane Results 
3 . . 1 Visual Features 

Visual features of the flame are the same as reported by 
Martin[38] and Yetter[i|23. Photographs of the combustor 
operating under the four different test conditions in the 
present study are presented in Figures 3-3 to 3-6. The visi- 
ble flame is small in diameter but very long, extending even 
past the sampling station, for the high co-swirl (55°) case. 
The flame diameter remains almost constant along its length 
and hence appears as a long, cylindrical, luminous core in 
the combustor. For the moderate co-swirl case (30°) the 
visible flame has similar features as for the high co-swirl 
case except that for 30° co-swirl the flame is shorter and 
the diameter of the flame near the entrance to the combustor 
is larger and gets smaller downstream. The flame looks like 
an extended bubble. In both the co-swirl cases the flames 
are smooth and appear to be less tu''bulent than in the coun- 
ter-swirl cases. 

The counter-swirl cases are characterized by greater noise 
and appear more turbulent. The flame surface is not smooth 
but is broken and shows rapid fluctuations indicating higher 
levels of turbulence. The flames are much larger in diame- 
ter, shorter in length and appear as luminous bubbles. The 
luminous bubbles are indicative of the recirculation zones. 
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3.^.2 Temperature arid Velocity 

The temperature and velocity profiles for the four test 
conditions listed in Table 2,3 are presented in Figures 3-7 
through 3-10, 

Profiles for the co-swirl oases (Figures 3-7 and 3-8) show 
steep temperature and velocity gradients. The high co-swirl 
case has a slightly higher centerline temperature and shows 
steeper gradients than the 30° co-swirl case. The axial vel- 
ocity profile for the 55° co-swirl case (Figure 3-7) shows 
almost a step change in axial velocity in the region of 
r/R^sO,5 indicating a hot and fast moving core of fluid con- 
fined to the center of the combustor and bears out the visual 
observation of a long cylindrical core of hot gases. The 
profiles for the 30 ° co-swirl case show gradients which are 
spread over a slightly larger radial thickness and are less 
steep (than for the 55° co-swirl case) indicating a greater 
mixing of the hot, high velocity central core with the outer 
annular jet. In each case the axial velocity gradient is 
confined to a narrow annular region (the shear layer). This 
shear layer is caused by the accelaration of the inner core 
due to reaction. The fact that there is such a shear layer 
which is confined to a narrow sheet in each of the co-swirl 
cases indicates that turbulent mixing across the combustor is 
poor. Turbulent mixing would smooth out the steep velocity 
and temperature gradients. The tangential velocity profiles 
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show high tangential velocities for the 55° case and moderate 
velocities for the 30° case as expected. 

The counter-swirl oases show more gradual radial varia- 
tions in temperature and velocities (Figures 3-9 and 3-10). 
Gradients are spread out across the entire radius indicating 
a high degree of trubulent mixing of the outer air with hob 
combustion product's. As a result, the temperatures near the 
centerline are much lower than in the co-swirl cases and the 
temperatures near the wall are higher than in the co-swirl 
cases. The axial velocities vary gradually across the whole 
radius of the combustor. A pronounced shear layer is not 
seen in the counter-swirl cases. In both the counter-swirl 
cases the outer swirl completely overpowers the inner swirl 
and hence the tangential veocities near the centerline are in 
the same direction as the tangential velocities in the outer 
flow. Lee*s[^^] cold flow measurements also demonstrate the 
dominance of outer swirl by the time the flow approaches 2 
diameters downstream of the inlet. 

It is seen from the temperature and velocity profiles (as 
well as from the composition profiles discussed in the next 
sub-section) that in- going from co-swirl to counter-swirl, 
turbulent mixing is enhanced. Two mechanisms may account for 
this increase: first, the tangential slip velocities across 
the inter-jet shear layer are increased in going from co- 
swirl to counter-swirl increasing turbulence generation in 
the shear layer; and second, swirl in the presence of a 
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strong positive radial density gradient, as in the co-swirl 
cases, suppresses turbulence production which in turn dampens 
the turbulent exchange of mass and momentum between the two 
jets . 

A comparison of Figures 3-9 and 3-10 shows that the axial 
velocity and temperature profiles are flatter and the center- 
line temperature is lower in the 30° counter-swirl case than 
in the 55° counter-swirl case indicating a greater mixing in 
the former case. This fact may be attributed to the greater 
suppression of trubulent mixing in the presence of density 
gradients due to the larger swirl in the high counter-swirl 
case . 

3.^.3 Exhaust Gas Composition 

Composition profiles on a dry basis are presented in 
Figures 3-11 through 3-1^. Composition results are also 
presented on a wet basis in Tables 3.2 through 3.5 along with 
the mass flux weighted average values for each case. The 
average values for all the oases are listed in Table 3.6 for 
comparison. Figures 3-11 through 3-1^ also show the profiles 
of the local equivalence ratio treated as a conserved 
scalar. 4 ^ is calculated from the local ratio of the total 
carbon atoms in the product species to the total oxygen atoms 
in the product species relative to the stoichiometric carbon 
to oxygen ratio. At the inletcj3=0.8 in the inner jet and 4^=0 
in the outer jet. Therefore, the profile of 4 ^ at the exit 
gives an indication of the extent of mixing between the inner 



and outer jots. The values of the overall equivalence ratio 

cb are also shown for each case, 
oa 

A significant feature of the exhaust gas composition data 
is that no intemediate lower hydrocarbons are present in 
detectable quantities. Lower hydrocarbons and alcohols are 
typically present in substantial amounts in the low tempera- 
ture oxidation or pyrolysis of hydrocarbons[79 , 80, 81 3 . 
Results of the present study, however, indicate that once the 
propane molcule is broken to initiate reaction, the reaction 
proceeds quickly to form CO or COg. Reaction quenching rates 
in the combustor are not high enough to freeze the reactions 
of intermediate hydrocarbon species. 

Composition profiles also show turbulent mixing trends 
similar to those shown by the temperature and velocity pro- 
files. Inspection of the 4* profiles (Figures 3~11 and 3-12) 
shows that the value of 4 is close to 0.8 near the center and 
falls steeply near r/R^sQ.^I in the co-swirl cases showing 
that there is little mixing between the inner and outer jets. 
For the 55° co-swirl case the temperature and equivalence 
ratio fall steeply near r/RQ=0.3 and 0.7. This region 
defines the extent of mixing between the inner and outer jets 
and will be called the mixing layer. The mixing layer is 
thicker for the 30° co-swirl case and extends between 
>^/Rq = 0.2 and 0.7* The mixing layers are denoted by arrows in 
Figures 3-11 and 3-12, The 4^ profiles for counter-swirl 
cases, Figures 3-13 and 3-1^, are much flatter and the <P 
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values are close bo the value across the entire radius of 
the combustor. Hence the mixing layer extends from the cen- 
terline to almost the wall of the combustor in the counter- 
swirl cases. The <j> profiles also show a slightly better mix- 
ing in the 30° counter-swirl case compared to the 55® 
counter-swirl case. 

The co-swirl composition profiles (Figures 3-11 and 3-12) 
have steep concentration gradients with little evidence of 
mixing with the high co-swirl case showing the steepest gra- 
dients. These gradients occur within the mixing layer bet- 
ween the hot inner jet of combustion products and the cold 
outer air jet as is to be expected. There is little or no 
fuel in the central hot core indicating combustion. Corres- 
pondingly, low concentrations of oxygen and large amounts of 
combustion products namely CO and COg are seen in the central 

core. The concentrations of CO and CO^ fall steeply while 

the concentrations of fuel and oxygen increase steeply across 
the mixing layer. CO concentration decreases less rapidly 
than the COg concentration in the mixing layer due bo the 

f 

quenching of CO oxidation reactions in the mixing layer as a 
result of dilution, as demonstrated by Fenimore and 
Moored?]. Outside the mixing layer, negligible amounts of 
CO and COg are found. The unburnt fuel in the exhaust comes 
from the outer edges of the inner inlet jet (mechanism is 

discussed in Section 3.8). Hence the maximum unburnt fuel 
concentration is expected in the mixing layer between the 

inner and outer jets. As a result of poor mixing in the co- 
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swirl cases, most of the unburnt fuel is confined within the 
narrow mixing layer. The steep decline in fuel concentration 
near the wall (r/R^sO.?) indicates the edge of the mixing 
layer. The steep decrease in fuel concentration towards the 
centerline is a result of combustion. 

The counter-swirl composition profiles are more uniform, 
with smaller gradients, and they show evidence of a greater 
degree of mixing. Substantial amounts of oxygen and unburnt 
fuel can be seen on the centerline as a result of greater 
mixing. The unburnt fuel in the counter-swirl cases also 
comes from the outer edges of the inner inlet jet, but it is 
distributed across the entire diameter at the exit as a 
result of mixing in the post combustion zone. The mixing, 
hov/ever, is not sufficient to produce very uniform concentra- 
tion profiles across the entire diameter of the combustor. 
The concentration of CO and of COg decrease towards the wall 
but some amount of CO and COg can be found close to the wall 
unlike the co-swirl cases. Fuel concentration decreases near 
the wall but not as steeply as in the co-swirl cases. 

The composition profiles also point to a slightly greater 
mixing in the moderate counter-swirl case than in the high 
counter-swirl case; oxygen and propane concentrations are 
higher on the centerline and their profils are much flatter 
for the moderate counter-swirl case. 
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3.^JI NO-NO^ 

Radial distributions of NO, NOj^ and NOg for the four test 
conditions are presented in Figure 3-15 through 3-18 as well 
as listed in Tables 3*2 through 3.5. The overall NO and NOj^ 
values (Table 3.6) ore quite low in agreement with previous 
results[26, 33, ^<0 , ^12] for lean premixed combustion. 

The NOj^ values follow the temperture values closely. Max- 
imum local values of NO^ occur in the co-swirl cases for 
which the maximum temperatures are observed, with high cen- 
terline NOj^ values decreasing rapidly to nearly zero in the 
mixing layer where temperatures are lower. NO values on the 

A 

centerline in the counter-swirl cases are nearly 10 times 
lower than in co-swirl cases due to the lower centerline 
temperatures in counter-swirl as well as due to the increased 
mixing in counter-swirl cases which causes the NO to be dis- 

A 

tributed across the entire diameter of the combustor. 

The importance of quenching in the NO to NOg conversion 
(discussed in Section 1.2,3) is clearly demonstrated by the 
profiles. In the co-swirl cases NO is predominantly NO in 

A 

the central core where quenching is minimal. Maximum amounts 
of NO are found in this region due to the high temperatures 
and reduced quenching. NO 2 , however, predominates in the 
mixing layer. The NO^ profiles peak in the mixing layer, 
with the peak being more pronounced in the 55° co-swirl case 
indicating a thinner mixing layer, than in the 30° case. 
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The NOj^ in the eountor-swirl oases is predominantly NOg 
over the entire radius at the exit. The large fraetions of 
NOg/NOjj can be attributed to quenching across the combustor 
in the post-cornbustion zone (see Section 1,2.3) by the mixing 
of the hot eombustion products and unburnt fuel with the 
outer air. 

Since large HOg concentrations are required for large 
NOg/NOjj fract^tons (Equation 1.2), the reduced levels of oxy- 
gen, which may limit HOg formation (Equation 1.3), may also 
be a factor in the reduced NOg values on the centerline in 
the co-swirl cases. 

Table 3.6 shows that inspite of locally high values of NO 
in the co-swirl cases the mass flux weighted average of NO^ 
values in the exhaust is of the same order for all the swirl 
conditions w.^ th slightly higher values for co-swirl cases. 
But, a larger fraction of NO^^ is NOg in the counter-swirl 
cases . 

It is believed that the use of cooled probes for sample 
withdrawl especially near reaction zones can result in the 
measurement of erronously high values for NOg concentration 
due to the reaction of NO with radicals in the sample probe 
druing the quenching process[30]. In the present study it 
can be concluded that the large fractions of NUg/NO^^ measured 
in the counter-swirl cases are real and not a result of 
quenching in the probe since large NOg/NO^^ fractions are not 
measured in co-swirl cases under the same probe quenching 
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conditions* Moreover, probe formed NOg is expected to be low 
in the present study since sampling is done well downstream 
of the combustion zone at a location where radical concentra- 
tions are expected to be low. 

3.^.5 Efficiencies 

Figure 3-19 shows the efficiencies obtained under the var- 
ious swirl conditions. The chemical efficiencies plotted in 
the figure are based on unburnt fuel with correction for CO 
In the exhaust. Efficiencies are also listed in Tabls 3.2 
through 3*6. 

Efficiencies are quite low since a lot of fuel is trapped 
in the mixing layer and in the outer flow and remains unburnt 
due to dilution and the resultant quenching. An important 
observation is that the efficiencies do not change signifi- 
cantly in going from coun ter-swirl to co-swirl though discus- 
sion in the preceding subsections suggested decreased mixing 
and therefore decreased quenching in the co-swirl cases. The 
almost constant efficiency seen in the present study is in 
contradiction to the trend observed by Yetter[^2] where the 
efficiency increased in going from counter to co-swirl . 
Hence additional results are required to explain the constant 
efficiencies observed in the present study. The additional 
experiments and their results are discussed in the following 
sections . 


It can be seen from Figure 3-19 that the thermal effi- 
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ciency is much Xower than the chemical efficiency for each 
case. Since the thermal efficiency is calculated based on 
the sensible enthalpy flux at the exit, the heat loss from 
the uninsulated combustor test section can account for the 
difference in chemical and thermal efficiencies. It is pos- 
sible to estimate the heat loss required to account for the 
difference as 

= m^q'j, ( \ -- \ ) (3.1^) 

It is seen that a heat loss of about 5 Kcal/s in each case 
can completely account for the difference in the efficien- 
cies. 

Heat can be lost from the hot gases by convection and gas 
radiation to the test section wall. This heat is transported 
by conduction across the wall thickness and is lost to the 
ambient. For the co-swirl cases the temperature profiles, 
Figures 3-7 and 3-8, suggest that convection to the wall is 
negligible (temperature profiles are flat near the wall; in 
fact the wall seems to be slightly hotter than the gas near 
the wall). However, due to the long luminous flame and high 
temperatures in the core, heat loss through radiation to the 
wall can be significant. In the counter-swirl cases the 
flames are shorter and the centerline temperatures are not 
very high near the exit. Hence radiation is important only 
near the inlet. The temperature profiles, Figures 3~9 and 
3-10, show gradients near the walls suggesting that convec- 
tive heat loss is important in the counter-swirl cases. Sira- 
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plified heat transfer analyses are done to estimate the mag- 
nitudes of the convective and radiative heat losses to ascer- 
tain whether a heat loss of 5 Kcal/s from the combustor is 
possible for each of the swirl conditions. 

Co-swirl ; Heat loss is almost entirely due to radiation. 
Combustion products, especially COg and water vapor, emit at 
infrared wavelengths[82, 833. Visible and ultraviolet band 
spectra of flames arise from electronic transitions of free 
radicals like OH, CH, Cg, CN and HCO in hydrocarbon flames. 
The reaction zone also emits continuous spectra, produced by 
recombination of dissociated molecules. For a first approxi- 
mation in calculations, the flame is often considered as a 
black body, and the fact that this is not entirely true is 
corrected for by an empirical factor smaller than 1[8^13. 

The flame in the co-swirl cases is approximated by a cyl- 
inder of radius 2.5 cm (from Figures 3-7 and 3-8) and length 
61 cm (upto the measurement station) and an average tempera- 
ture of 1800°K. The wall temperature is 3^0°K. Since the 
flame is a clean hydrocarbon flame the emissivity will be 
considerably lower than 1 and hence a value of 0.3 for the 
emissivity is reasonable to assume. The Stefan-Bol tzman law 
is used to calculate the radiative loss, with the surface 
area obtained from the flame dimensions mentioned above. The 
heat loss is found to be 4.1 Kcal/s. It should be noted that 
the flame temperatures near the inlet are much higher and the 
heat loss can be greater than the 4.1 Kcal/s calculated 
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above. Therefore, it is possible to account for the heat 
loss of approximately 5 Kcal/s for the co-swirl cases. 

Counter-swirl ; Heat transfer is primarily through convec- 
tion. An estimate for the turbulent thermal conductivity of 
the gas near the wall, k’ , is obtained from k'= <?c u’l, with 
turbulent Prandtl number=1, where u’ and 1 are the turbulent 
velocity and length scales respectively. Measurements of 
mean velocities made by DepskyCYl] show that u’ is of the 
order of 3 m/s. 1 is taken to be 0.005 m which is approxi- 
mately the thickness of the boundary layer as see'.i from the 
mean axial velocity profiles for the different swi.. 1 condi- 
tions. The values of Cp and ? at 350°K are 0.24 Kcal/kg 
and 0.998 kg/m^ respec tiv ely [ 84 ] , which give a value of 
k'=3-6x10 ^ Kcal/s m °K. Heat loss by convection to the wall 
is given by Q^= k’A(dT/dr), where A is the inner surface area 
of the test section. (dT/dr) is estimated to be 7.4x10^ °K/m 
from Figure 3-9 and is found to be 5.22 Kcal/s. It should 
be noted that the temperature gradient at the wall near the 
inlet will be lower and the convective loss will be lower. 
But the radiation losses are higher near the inlet. Hence it 
is reasonable to expect a heat loss of about 5 Kcal/s in the 
counter-swirl cases also. 

As an additional check, the average temperature difference 
across the thickness of the quartz test section required to 
lose 5 Kcal/s of heat is found to be 28°C (Thermal conductiv- 
ity of quartz = 2.2x10“''^ Kcal/s m°K, 


wall thickness = 



73 


0.0025 m) . It is reasonable to expect this temperature 
difference in actuality since the outside of the test section 
is exposed to the ambient air (295°K) while the inner wall is 
exposed to a temperature of about 3^0®K at the exit. How- 
ever, actual measurements have not been made to determine the 
exact temperature difference across the wall thickness. (The 
ambient temperature is marked in the Figures 3-7 to 3-10 
along with the temperature profiles) . 

It should also be noted that the estimate of 5 Kcal/s heat 
loss to account for the difference in the chemical and ther- 
mal efficiencies is based on the assumption that the chemical 
efficiency is the correct efficiency. However, air leakage 
into the samples during composition measurements results in a 
lower value for the mean unburnt fuel fraction and hence the 
chemical efficiency are overestimated. Therefore, the actual 
heat loss required is less than 5 Kcal/s. It is also to be 
noted that the thermal efficiencies would be about 2 % higher 
than what are shown in Figure 3-19 if the temperature pro- 
files are corrected for radiation losses from the thermocou- 
ple. 

It is concluded that the difference between the calculated 
chemical and thermal efficiencies is due to the heat loss 
from the combustor. It is seen, based on the chemical effi- 
ciency, that the heat loss of 5 Kcal/s is about 25 % of the 
heat generated in the combustion process. It is possible 
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that this heat loss by itself could partly account for the 
low efficiencies observed in the present study. 

3 • 5 Comparison with Previous Results for Methane in the 
Combustor 

Yetter ’ s[ and 0ven's[40] results with methane fuel in 
the combustor are compared in this section. 

Temperature and velocity results in the present study are 
qualitatively similar to Tetter's results which also show 
steep gradients for co-swirl cases (Figures 3-^a and 3~5a in 
C^2]). The maximum (centerline) temperatures in the co-swirl 
cases in the present study, after the addition of radiation 
correction of approximately 90‘^K, are within 1.5^ of the max- 
imum temperatures observed by Yetter. The absolute values of 
the velocities in Tetter's study are much higher than in the 
present study (See Figure 3-26). Since the same probe was 
used in both the studies and the calibration done by Yetter 
was also in a jet of low turbulent intensity, his velocity 
data need to be corrected as described in Section 2.3* 

Tetter's results (Figures 3-1a and 3-2a in [42]) also show 
a slightly lower centerline temperature and flatter- tempera- 
ture and velocity profiles for 30° counter-swirl case than 
the 55° counter-swirl case indicating better mixing in the 
30° case. The 55° counter-swirl centerline temperature in 
the present study is significantly higher than in Tetter's 
results, but the wall temperatures in the present study espe- 
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daily in the oounter-swirl cases are lower* 

A significant observation is that the maximum temperature 
and velocity gradients in Yetter's results are found at lar- 
ger radial distances than in the present study at correspond- 
ing axial stations, indicating that the central hot core is 
larger in diameter in Yetter’s tests. 

Though a number of similarities in composition profiles 
can be seen between Yetter’s and the present results, there 
are some significant differences, Yetter’s composition 
results show that the fuel concentrations do not peak near 
r/R^=0.5, rather they peak at the combustor wall. This fact 
seems bo imply a significant difference between the way the 
fuel is mixed and burnt under his test conditions and the 
present test conditions. Yetter's results also show pro- 
nounced CO peaks in the mixing layer indicating quenching of 
the CO oxidation reactions in the mixing layer. The relative 
locations of CO gradients in the two studies also indicates 
that the mixing layers in Yetter's study are at larger radial 
distances and hence the central hot core is larger in diame- 
ter . 

There are no significant differences in the NO-NO^ results 
between the present study and Oven's results. Oven's results 
for high co-swirl (Figure 16 in [40]) also show low NO /NO 

^ X 

ratios in the central core and an NO^ peak in the mixing 
layer while almost all of the NO is NO in the counter-swirl 

X 

case (Figure 18 in [40]). Average NO^ values in the exhaust 
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are not available from previous resultsC^O, ^2] for compari- 
son. Though 0ven[^0] has not made velocity measurements, the 
radial measurements of temperature and composition at various 
axial locations upto 6 D^ downstream from the combustor inlet 
show trends which are consistent with observations at 6 D 

o 

made in this study. 

The largest differences between results of the present 
study and Tetter’s results are found in the efficiency 
values. Firstly, the efficiencies in Tetter’s study are, in 
general, higher and secondly they increase significantly in 
going from counter-swirl to co-swirl (Figure 3-19). This 
trend is a serious deviation from the nearly contant effi- 
ciencies seen in the present study. The only significant 
differences in the operating conditions between the two stu- 
dies are in the inner swirl number and axial velocity ratio. 
Tetter used Sj^ = 0.523 and U^/U^=1.5 for his study. The 
slightly lower axial velocity ratio (= 1 . 3 ) in the present 
study should, if anything, increase the efficiency according 
to the study of the effect of axial velocity ratio conducted 
by Tetter[il2]. Effects of changes in inner swirl on combus- 
tor performance have not been studied so far. 

The chemical and thermal efficiencies in Tetter’s study 
are in good agreement with each other indicating a lower heat 
loss than in the present study. It is possible that the 
erronous high values for the velocities reported in Tetter’s 
study, which contribute to an increase in the sensible 
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enthalpy flux, have caused overestimates of the thermal 
efficiencies, thereby underestimating the heat loss. 

3.6 Phase 2 ; Cold Flow Measurements 

Differences in the fuel distribution in the exhaust plane 
between Tetter's results and the present study suggest dif- 
ferences in the way the fuel is mixed and burnt in the com- 
bustion zone. 0ven[^}0] has performed measurements at the 
inlet under cold flow conditions to test the efficacy of 
fuel-air premixing in the inner jet. His results indicated 
that the maximum fuel concentration occurs on the centerline 
and concentrations at the edge of the inner jet are 3555 lower 
than the maximum (Figure 28 in [^0]). One of the causes for 
such a distribution was thought bo be the centrifugal forces 
in the swirling jet differentially forcing the air, which has 
a higher molecular weight than methane, towards the outer 
edge of the jet. Such a distribution was considered benefi- 
cial to the combustor efficiency since less methane would be 
diluted at the outer fringes of the inner jet in the inlet 
region of the combustor. 

As a first step in the cold flow measurements, in the pre- 
sent study, the fuel distribution at the exit was measured 
for one of the operating conditions. Figure 3-20 shows the 
radial distribution of propane at the exit plane (6D^) in 
cold flow under the 30° counter-swirl condition. The figure 
shows that the extent of mixing is not sufficient to produce 
a uniform concentration of propane across the exit. Propane 
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concentration at the wall is 505C lower than the maximum con- 
centration occuring on the centerline. The profile rein- 
forces the observations in the hot flow of decreasing fuel 
concentrations near the combustor wall. 

Possible differences in the injection mechanism and pre- 
mixing in the inner jet due to different molecular weights of 
propane and methane were examined by inlet plane measure- 
ments, Results of these measurements, presented in Figure 
3-21, show that the mixing of propane and methane in the 
inner jet are not significantly different. The inlet distri- 
bution of the fuels in slightly different from that observed 
by Oven in that the concentration at the edge of the inner 
jet is only 20 % lower than on the centerline in the present 
study, in comparison to 35 % in Oven’s study. This may mean 
that more fuel could be diluted at the outer fringes of the 
inner jet causing lower efficiencies in the present study. 
But the differences in the variation of efficiency with outer 
swirl still remains to be explained. 

The cold flow measurements suggest that the differences 
between the efficiencies observed by Yetter and in the pre- 
sent study are not likely to be due to differences in the 
premixing of the fuels propane and methane. Measurements 
were therefore performed with methane fuel in the combuslor 
for similar test conditions as for propane in order to -^^om- 
pare directly with Tetter’s results and to determine if the 
differences in the efficiencies between the two studies are 
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due to the differences in the chemistry of combustion of 
propane and methane. 

Phase 2 ; Methane Results 

Results are presented for the two test conditions listed 
in Table 2.3. Temperature and velocity profiles are pre- 
sented in Figures. 3~22 and 3-23. Composition profiles are 
presented in Figures 3-24 and 3-25. Concentrations in ppm on 
a wet basis, the average values and the efficiencies are 
listed in Tables 3.7 and 3.8. Average values are compared in 
Table 3.6. NO-NOj^ measurements were not performed for meth- 
ane . 

The temperature and velocity profiles are almost identical 
to the profiles for propane for corresponding swirl condi- 
tions (Figures 3-7 and 3-9). The maximum temperature and 
velocity gradients occur at identical radial locations. 

Composition profiles for methane are qualitatively identi- 
cal to the profiles for propane. Steep concentration gradi- 
ents (especially CO gradients) occur at the same radial loca- 
tions for propane and methane for the respective swirl 
conditions. Hence all the remarks in the discussion of pro- 
pane results in Section 3.4 are pertinent to methane results 
also, with the exception that the CO peaks are more pro- 
nounced in the methane results. 

Efficiencies for the methane case (Figure 3-19) are 
approximately 5% higher, but as in the case of propane the 
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efficiencies do not change significantly in going from coun- 
ter to co-swirl. The temperature profile corrected for radi- 
ation losses from the thermocouple for the high co-swirl 
methane case is shown and compared with the uncorrected temp- 
erature profile in Figure 3-22. Only the high temperatures 
close to the centerline ore altered by the radiation correc- 
tion. The value of the thermal efficiency goes up by only 2 % 
when the corrected temperatures are used for the calculation. 
This increase no doubt brings the value of thermal efficiency 
closer to the respective chemical efficiency, but does not 
alter the conclusion that efficiency is not significantly 
changed by variation of outer swirl for the test conditions 
studied . 

A comparison of temperature and velocity profiles reported 
by Yetter and in the present study is shown in Figure 3-26. 
It is apparent from Figure 3-26 that for Tetter's operating 
condition, the central core is hotter and larger in diameter 
which strengthens the same observation made in Section 3*5. 

The value of inner swirl number for Yetter (S^=0,523) is 
only slightly different from S^ = o,il95 used in this study. 
However, even a small change in seems to affect the 
combustor operation substantially. 

3. 8 Proposed Mechanism for the Combustor Operation 

The following mechanism for the combustion operation is 
proposed which will consistently explain the observations in 
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the present study as well as Tetter’s study, 

It is proposed that most of the fuel burnt is burnt in 
front and on the sides of the recirculation zone, in a region 
very close to the combustor inlet. Fuel from the inner 

regions of the inner jet is burnt in this region of low to 
moderate velocities, thus forming the intense reaction zona 
(shown in Figure' 1-1) which will be referred to as the 
'flame’ in the ensuing discussion. (It has been estab- 

lished[i}7 , 71 ] that the role of the recirculation zone is to 
stabilize the flame by providing a region of low velocity 
flow in front of the recirculation zone) . Inefficiency in 
the combustor is caused by the .inability to burn the fuel 
flowing at the edges of the inner jet which flows around and 
at large radial distances from the recirculation zone and is 
lost to the combustion process. High efficiencies result 
from the ability of the flame to propagate radially outward 
in order to burn the fuel flowing in the outer edges of the 
inner jet. 

The following four characteristics are considered impor- 
tant in determining how much fuel is burnt in front of the 
recirculation zone and the ability of the flame bo propagate 
radially: 1) The size of the recirculation zone, 2) inlet 

turbulence, 3) velocities around the sides of the recircula- 
tion zone and ^) the local equivalence ratio on the sides of 
the recirculation zone. 
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Boyler's[47] mean CH radiant emission data are plotted, 
with the reverse flow region as measured by Depsky[713» in 
Figure 3-27. The elliptical arc represents the reverse flow 
region (the other symbols in the Figure are explained later 
in this section). The extent of the recirculation zone is 
about the same as that of the region of reversed flowC713. 
It is seen from Figure 3-27 that the recirculation zone is 
larger for the counter-swirl case than for the co-swirl case. 
Since a larger recirculation zone has a larger region of 
reduced velocity in front of it, it is proposed that a larger 
amount of fuel is burnt in front of the recirculation zone in 
the counter-swirl case than in the oo-swirl case. The above 
hypothesis is necessary to explain the constant efficiencies 
observed in the present study as will be seen later in this 
section. A rigorous check on this hypothesis would be to 
clearly delineate the reaction zone in front of the recircu- 
lation zone, draw streamtubes from the inlet jet and deter- 
mine what mass of fuel-air mixture enters the reaction zone 
for each of the outer swirl conditions. Unfortunately, suf- 
ficient data near the inlet is not available at present to 
make this check. 

A recent study by Moreau and Borghi[85] of a turbulent 
premixed flame shows large changes in combustion efficiency 
for small changes in inlet turbulence. Lee’s[44] measure- 
ments in the present combustor for cold flow show large r.m.s 
values for axial and tangential velocities near the inlet for 
counter-swirl compared to the co-swirl case, which might also 
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contribute to a larger amount of fuel being burnt in front of 
the recirculation zone in the counter-swirl case than in the 
co-swirl case. It should be noted that in Lee’s[M43 cold 
flow measurements there is no recirculation zone for the co- 
swirl case, whereas under combusting conditions there is a 
recirculation zone for the co-swirl case also. Depsky’s[71] 
measurements under combusting conditions show large values of 
r.m.s axial and tangential velocities at the inlet for co- 
swirl also. It is worthwhile to make more detailed measure- 
ments of the inlet turbulence (in front of the recirculation 
zone) to characterize its effect on the flame. 

The ability of the flame to propagate radially outward 
depends greatly on the local velocity and equivalence ratio 
on the sides of the recirculation zone. In Figure 3-2? the 
circles represent the peaks in the mean CH radiant emission 
profiles and the squares represent the half maximum points 
from Reference 47. The triangles represent the location of 
the axial velocity peaks and the vertical bars represent the 
region of steep radial gradients of axial velocity 'from 
Reference 71 for nearly the same operating conditions includ- 
ing the same inner swirl number as in Reference 47 and in the 
present study. The peaks in the axial velocity (Figure 3-27) 
are, in part, a result of the acceleration of the flow around 
the recirculation zone as opposed to combustion induced 
acceleration and are seen in Lee’ s[44]‘ cold flow measurements 
as well. Chemiluminescent emission from CH has been estab- 
lished as a signature of the reaction zone in the combus- 
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tort 07 3. Figure 3-S7 shows that the maxima of meau GH emis- 
sions oQour outside the reoirouXatton zone but well interior 
to the region of peak axial velocity. This implies that the 
flame is not able to propagate aoross tne region of high vel- 
ooities. The peaks in the axial veXooities are as high as 
1.25Uj^ where is the average axial velocity in the inner 
jet. The flame is able to propagate radially until the steep 
gradients in velocity are encountered. Then the radial pro- 
pagation is slower since the flame is almost parallel to the 
flow direction due to the high flow velocities. This slower 
radial propagation as the flow moves downstream is possible 
only if the loeal equivalence ratios are above the local 
inflammability limits determined by the local turbulence, 
temperature and velocity conditions. For both the oo- and 
counter-swirl oases the peak velocities are nearly the same 
(1.25U^) near the inlet but the peak value is slightly lower 
and the peak is broader by the time the flow reaches 1D^ 
downstream as a result of increased momentum transport due to 
greater mixing in the ooun ter-swirl case (Figures 7 and 11 in 
[713). In the oo-swlrl oases, however, the local equivalenoe 
ratios in the region of the steep gradients of axial velocity 
are expected to be higher than the local inflammability lim- 
its due bo reduced mixing with the outer jet. Therefore, the 
flame can propagate radially outward though it does so very 
gradually on aooount of the high axial velocities and the 
flame being almost parallel to the flow. In oounber-swirl , 
in addition bo the velocities on the sides of the reciroula- 


tion zone being high, the looal equivalenoe ratios in that 
region are expected to decrease rapidly as we go downstream 
due to increased mixing with the outer air. Hence the flame 
cannot propagate radially outward as the flow moves down- 
stream. Figure 3-27 shows that the flame as denoted by the 
peaks of mean CH emissions does not propagate radially out- 
ward to the excisnt it does in the oo-swirl case. The greater 
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the radial propagation of the flame, the larger the number of 
streamtubes it is likely to out, thereby burning more fuel. 

Based on the mechanism discussed above, a smaller amount 
of fuel is burnt in front of the smaller recirculation zone 
in the oo-swirl case than in the counter-swirl case. But more 
fuel is burnt on the sides of the recirculation zone in the 
co-swirl case than in the counter-swirl case. 

In the present study, it so happens that the amounts of 
fuel burnt in front and on the sides of the recirculation 
zone add up to approximately the same amount in the co- and 
counter-swirl cases resulting in the same efficiency. As 
seen from the discussions above, this insensitivity of the 
efficiency to outer swirl does not imply that the combustion 
mechanism is unaffected by outer swirl. It could be consid- 
ered fortuitious that the operating conditions set in the 
present study resulted in a constant efficiency with varying 
outer swirl. On the other hand, with flow conditions result- 
ing in a smaller recirculation j:one in co-swirl than in the 
present study, one can possibly achieve higher efficiencies 
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for counter-swirl than for co-swirl though there is more 
dilution from the outer air in the counter-swirl than in the 
co-swirl case. There are also operating conditions for which 
the efficiencies for co-swirl are greater than those for 
counter-swirl. 

Yetter’s[42] observation of a higher efficiency for 
co-swirl than for counter-swirl may be explained as follows. 
A study of the size of the recirculation zone formed from a 
single jet for varying swirl numbers made by Syred and 
Beer[2] shows that the recirculation zone increases in size 
with increasing swirl numbers. This suggests that for 
operating conditions in [4^3, the recirculation zones were 
probably bigger than in the present study due to the larger 
inner swirl in [42]. This means that in co-swirl cases in 
[42], a larger amount of fuel is burnt in front of the 
recirculation zone as well as that the flame propagates 
farther radially than in the present study since the high 
velocity peaks occur at larger radial distances due to the 
larger recirculation zone resulting in much higher 
efficiencies (Figure 3-19) and a larger central hot core 
(Figure 3-26) than in the present study. The increase in 
efficiency in the counter-swirl cases is not much (Figure 
3-19). There seems to be an upper limit to the efficiency in 
the counter-swirl cases as a result of dilution from the 
outer air. Beyond a certain size, an increase in the size of 
the recirculation zone in the counter-swirl cases is of no 
use, since, eventhough the flow conditions are suitable for 
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more oombusfcion> the Xonal equivalence ratios are too low as 
a result of dilution. 

Therefore, it appears that even small changes in inner 
swirl can cause large changes in efficiency. Exploratory 
measurements performed by the author, the results of which 
are not presented here for want of systematic data, show that 
minor adjustments of the butterfly valve controlling the 
inner swirl causes changes in centerline temperature at the 
exit of the order of 150°K. The experiment done for the mod- 
erate co-swirl case demonstrates that the combustor is very 
sensitive to inner swirl. 

The oxidation of CO and NO in the combustor to COg and NO^ 
respectively is explained as follows. The CO to COg and NO 
to NOg conversions being slower reactions are determined by 
conditions downstream of the combustion zone in the post-com- 
bustion zone. Inspection of the time mean isotherms in the 
combustor measured by Oven, et al.Cmj (Figure 3-28) shows 
that cooling rates on the centerline for the counter-swirl 
case can be as high as 250°K per millisecond (an axial veloc- 
ity of 50 m/s is assumed) . This cooling rate is much higher 
than the minimum cooling rate of 70*^K per millisecond at 
1270°K specified by Fenimore and Moore[17] for CO remaining 
permanently quenched. There is little or no cooling of the 
hot core in the co-swirl case (Figure 3-28). Hence large 
amounts of CO are found at the exit in the counter-swirl 
cases compared to the oo-swirl cases. Inspection of the time 
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mean isoplefchs of CO ineasuned by Oven, eb al.[M13 and 
presented in Figure 3-29 shows CO being rapidly converted to 
COg in the central core of the post-combustion zone for the 
oo-swirl case whereas large amounts of CO are quenched in 
this zone and mixed across the combustor for the counter- 
swirl cases. Large fractions of counter-swirl 

cases are a result of the increased quenching in the post- 
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combustion zone in the counter-swirl cases promoting the for- 
mation of H0« which oxidizes NO to N0„. 


CHAPTER 


CONCLUSIONS 

The nieasuretnents reported in this work in combination with 
results from previous measurements in the swirl combustor 
lead to the following major conclusions and suggest future 
experiments designed to understand specific aspects of the 
combustor operation. 

The exhaust from the combustion of propane in the combu- 
stor is free from lower hydrocarbons like ethane, ethylene 
and methane. 

There is no difference in the combustor operation for the 
fuels propane and methane. The use of either of the fuels 
results in nearly the same temperature and velocity profiles 
and the same efficiency for a given operating condition sug- 
gesting that differences in the chemistry of combustion bet- 
ween propane and methane do not affect the overall combustor 
operation. 

Practically all the combustion occurs in the region close 
to the combustor inlet upto about 1 * 50 ^ from the inlet, 
called the combustion zone. The CO to COg and NO to NOg con- 
versions occur in a region downstream of the combustion zone 
called the post-combustion zone. 

Turbulent mixing in the combustor is enhanced in going 
from high co-swirl to counter-swirl with the 30 ° counter- 
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swirl showing the greatest degree of mixing. The mixing 
trends are a result of interactions between tangential shear 
between inner and outer jets and swirl in the presence of 
positive radial densit;;?' gradients. 

The NOj^ levels in the exhaust are low as a result of lean 
preinixed combustion. Large fractions of NOg/NOj^ are found as 
a result of quenching in the post-combustion zone due to 
dilution from the outer air jet especially in the counter- 
swirl cases. Quenching of CO oxidation reactions in the 
post-combustion zone is the cause for the large amounts of CO 
in the counter-swirl cases. 

Combustion efficiencies in the present study are low and 
remain constant with varying outer swirl. A mechanism for 
the combustor operation is proposed which explains the cons- 
tant efficiencies in the present study as well as the 
increase in efficiency in going from counter to co-swirl 
observed by Yetter[i|2]. The flame (intense reaction zone) is 
stabilized in front of the recirculation zone. Inefficiency 
in the combustor stems from the inability of the flame to 
propagate radially outward and burn the fuel from the outer 
edges of the inner jet. The following four fluid-mechanical 
characteristics are considered important in determining the 
efficiency of the combustor: 

1) The size of the recirculation zone 

2) Inlet turbulence 

3) Velocities around the recirculation zone 
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^0 Local equivalence ratio on the sides of the 
re lirculation zone. 

It is apparent that these flow characteristics should be 
investigated in detail and related to the operating 
parameters of the combustor as a viable means to control 
combustion . 

It is evident that the combustor operation is very sensi- 
tive to even small changes in the inner swirl number. It is 
therefore necessary to study the effect which the inner swirl 
has on the above mentioned four flow characteristics for 
different outer swirl conditions. 

A few measurements are suggested for the future aimed at 
achieving the objectives mentioned above. It is necessary to 
delineate the flame precisely and draw streamtubes for the 
flow in the combustion zone to determine how much of fuel is 
burnt where in the combustion zone. It is important to use 
non-intrusive optical techniques for the measurements. Laser 
Doppler Velocimetry is suggested for making extensive mean 
and r.m.s velocity measurements especially in front of the 
recirculation zone. Temperature measurements using the Ray- 
leigh scattering technique will provide information on the 
local densities. The local velocity and density information 
can be used to draw streamtubes based on mass flow. Measure- 
ments of CH emissions using Emission Spectroscopy can be used 
to delineate the flame. Measurement of major species concen- 
trations in the cambustion zone using techinques such as 
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Raman Spectroscopy provides information on the local 
equivalence ratios. These measurements should performed 
for varying invfer and outer swirl conditions. Effects of 
other operating conditions like inner jet equivalence ratio 
and inner to outer jet axial velocity ratio on the combustor 
performance can also be studied. There are, however, numer- 
ous problems in the practical utilization of the suggested 
techniques like signal processing, noise, biases in measure- 
ment, maintaining uniform scattering cross-sections for Ray- 
leigh scattering measurements, spatial resolution, signal to 
noise ratio especially for low concentrations, laser beam 
quality and reliability, etc. Most of these problems are 
being overcome and the techniques are becoming available for 
accurate and reliable measurements. Optical access to the 
region of interest in the combustor is necessary. Figures 
3-3 to 3-6 show that the combustion zones extend into the 
inner tube of the combustor posing the problem of optical 
access to that region. But in general, the configuration of 
the present combustor provides good optical access and is 
well suited for optical techniques. 

Measurements similar to those suggested above can be per- 
formed in the post-combustion zone to study the dynamics of 
the CO to CCg and NO to NOg conversions. Conventional sam- 
pling and velocity probes are nevertheless still effective 
tools for analysis in regions far away from reaction zones 
like near the exit. New techniques using acousto-optic fil- 
ters (SMART sensors) [86, 873 for absorption spectroscopy are 
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being developed which may prove very valuable in combustion 
research on account of their speed and ruggedness of opera- 
tion, The development and application of these techniques 
will help in understanding the combustor better, 

Some modifications to the coK‘^ustor design could be 
attempted, For example, a shield can be provided around the 
fuel injector so that fuel is injected only into the central 
region of the inner jet. This will improve the combustor 
efficiency since little or no fuel is diluted on the outer 
edges of the inner jet. Another modification could be to put 
a diffuser-like ,'f»deld over the inner jet and the recircula- 
tion zone to delay the mixing of the outer jet with the inner 
jet. With a better understanding of the present combustor, 
suitable designs for practical applications of premixed pre- 
vaporized swirl combustors can be developed. 
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Table 2»1« 

Compositions (in ^ by volume) of 
calibration gas mixtures 


Mixture No. 


^2 

CO 

COg 

CH 4 


(air) 
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79.1 


5.0 

72.38 

7.02 

15.6 
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C 3 H 8 


10.0 


1.0 


Response factors and estimated errors 
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Test Conditions 


Outer Swirl 
PROPANE 1 

^0 

Ui 

Cm/s] 

u./u 

r 0 

“oa 

[m/s] 

‘(’i 

^oa 

55 ®Co - 

+ 0.56 

30.0 

1.32 

24.3 

0.81 

0.214 

30°co - 

•♦■0.20 

30.3 

1.25 

25.6 

0.80 

0.210 

30 °Ct - 

- 0.20 

30.1 

1.26 

25.3 

0.80 

0.210 

55 °Ct - 

-0.56 

30.2 

1.37 

23.9 

0.80 

0.224 

METHANE 1 

55 °Co - 

■♦• 0 . 56 

31.5 

1.33 

25.4 

0.81 

0.218 

30 °Ct - 

- 0.20 

31.6 

1.28 

26 . 2 

0.80 

0.213 


= 0.495 = 3.4 SCFM (Propane) 

= 8.5 SCFM (Methane). 


Table 3 . 1 . 


)0H 


0 values on either side of the centerline 
for the 55° co-swirl, propane case 
showing axisyininetry at the measurement station 
and reproducibility of combustor operation 


r( cm) 

e(°) 

fir ) 

sepp-rate : 

5.0 

-19.5 


4.0 

-22.5 

-22.0 

3.0 

-21.0 

-20.0 

2.0 

-16.0 

-15.5 

1.0 

-10.5 

-10.5 

0.5 

-6.0 

-6.0 

0.0 

0,0 

0.0 

0.5 

6.0 

5.5 

1.0 

11.0 


2.0 

16.5 


3.0 

21.0 


4.0 

23.0 


5.0 

20.5 
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Composition (wet basis), velocity, temperature, their average valves 
and efficiencies for the 30 ° co-swirl, propane case 
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CHEMICAL EFFICIENCY (UNBURNT FUEL ONLY) = 51.773% 
CHEMICAL EFFICIENCY (UNBURNT FUEL AND CO)= 51.439% 


Composition (wet Pasis), velocity, temperature, their average valves 
and efficiencies for the 30 ° counter-swirl, propane case 
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CHEMICAL EFFICIENCY (UNBURNT FUEL ONLY) = 57. 131Z 
CHEMICAL EFFICIENCY (UNBURNT FUEL AND CO)= 55. 202X 


Composition (wet basis), velocity, temperature, their average valves 
and efficiencies for the 55° counter- swirl, propane case 
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CHEMICAL EFFICIENCY (UNBURNT FUEL ONLY) = 54. 002X 
CHEMICAL EFFICIENCY (UNBURNT FUEL AND CO)= 51. 419X 
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2— 1j Schfeina."fcic diagram iihe swirl confbusijor tes1: facility. 
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Figure 2-2i Variable swirl generators 
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Figure 2-3. Cylindrical prote for sampling and velocity measurements 
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Figure 2-4i Sample train and analysis system. 
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(A) 


»omple injection 


4 port volve 


— 

poropok 0 



VW^ 
corbosieve S 


to detector 


(B) 

4 port volve 



to detector 


I Schematic diagram of column arrangement and 
flow paths in the gas chromatograph for valve 
positions A and B. 




Figure 2-6. Sample chromatogram. (The numbers on the sides 
of the peaks indicate detector attenuation. 

. The sample injection point (I) and points of 
valve switching to positions A and B are in- 
dicated. The time axis is compressed between 
peaks. The integrator trace at the bottom 
shows the number of line crossings). 
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b) One Cl the pressure tops is disconnected 
to meosure P 59 - Poim 


flow direction 



c) The probe is rototed by 59® to meosure 

Pt - P * 

T otm 

The difference between (b) ond (c) gives 


Figure 2-7i Velocity measurement scheme. 
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Figure 2-9* Pressure distribution around the surface of the cylindrical 
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Figure 3-1 i Mean temperature profiles for the 55° co-swirl, 
propane case showing axisymmetry in the com- 
bustor. 
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Figure 3-2 i Mean temperature profiles for two separate 
runs of the 30 ° co-swirl case showing repro- 
ducibility of combustor operator. 
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Figure 3-3i Photogragh of the combustor in operation' 
under 55° co-swirl, propane firing. 



Figure 3-4i Photograph of the combustor in operation 
under 30 ° co-swirl, propane firing. 
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Figure 3 - 5 « Photogragh of the combustor in operation 
under 30® counter-swirl, propane firing. 



Figure 3-61 Photograph of the combustor in operation 
under 55^ counter-swirl, propane firing. 




ORIGINAL PAGE IS 
OF POOR QUALITY 


127 



r/Ro 


Figure 3~?t Mean temperature and velocity profiles for 
the 55° C0“swirl, propane case. 
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Figure 3-8 i Mean temperature and velocity profiles for the 
30 ® co-swirl, propane case. 
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□m = 52.90 m/s 
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Figure 3-9 « Mean temperature and velocity profiles for the 
30 ° counter swirl, propane case. 





0.5 


MiQINAL PAGE m 
OF POOP QUAUTY 


1144 

51.60 m/s 
16.96 m/s 




1.0 0.8 0.6 0.4 0.2 0.0 0.2 0.4 0.6 0.8 1.0 

f /Rq 

Pig^ire 3-10i Mean temperature and velocity profiles for the 
55° counter-swirl, propane case. 
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Figure 3-11 i Mean composition and equivalence ratio profiles 
for the 55^ co-swirl, propane case. 
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Figure 3-12 i Mean composition and equivalence ratio pro 
files for the 30 ° co-swirl, propane case. 
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Figure 3“13« Mean coraposition^and equivalence ratio pro- 
files for the 30 ° counter-swirl, propane case. 
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Figure 3-15 « Mean concentration profiles of oxides of 

nitrogen for the 55° co-swirl, propane case 
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Figure 3-I61 Mean concentration profiles of oxides of 

nitrogen for the 30° co-swirl, propane case 
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Figure 3-17* Mean concentration profiles of oxides of 
nitrogen for the 30° counter-swirl case. 
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Figure 3-19* Chemical and thermal efficiencies for the 

different outer swirl conditions under propane 
and methane firing in the present study com- 
pared with those in Yetter’s [42] study. 







Figure 3 - 20 i Mean concentration profile for propane for 
the 30^ counter-swirl case under cold flow 


conditions . 
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Figure 3-21 i Fuel distribution in the inner jet at the 
inlet. Molefractions plotted are non- 
dimensionalized by their maximum values 
in each case . 
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Figure 3-22 1 Mean temperature and velocity profiles for 

the 55° co-swirl, methane case. Temperature 
values after radiation correction are also 
shown. 
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Figure 3-23 i Mean temperature and velocity profiles for 
the 30 ° counter swirl, methane case. 
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Figure 3-24j Mean composition profiles for the 55 ° co-swirl 
methane case. 
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Figure 3“25» Mean composition profiles for the 30 
counter-swirl, methane case. 
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Figure 3"26i Mean temperature and axial velocity profiles 
for the 55° co-swirl, methane case in the 
present study compared with the corresponding 
profiles in letter's L^-2] study. 
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Figure 3-28 i Kean isotherms (°K) in the combustor (from Reference m) 




APPENDIX A 


Equablona for Determining Combustor Operating Conditions 


Subscripts ; 
i - inner jet 
o - outer jet 
oa - overall 
f - fuel 

1. Inner Jet Mass Flow Rate (Ibm/inin) 

</FP7T)^ 

Where « Constant 24.1285 (determined from pitot 

tube measurements) 


(A.l) 


h “ differential pressure across the annubar 

(inches of water column) 


p “ static pressure at the annubar (psia) 


T = temperature at the annubar (^R) 


2. Outer Jet Mass Flow Rate (Ibm/min) 

= A 2 (>/Hm)o <A.2) 

Where A 2 - constant = 59.5491 (determined from pitot 

tube measurements) 

3. Volumetric Air Flow Rate (CF>l) 

Q = rii/p (A. 3) 


150 


* 


151 

Where p * density in the corresponding Jet (Ibrn/ff^) 

4. Inner Jet Mean Axial Velocity (ft/sec) 

- (Qi -f Qf)/A^ 

- 0.8692 [8.916(»W)?)i Qfl (A. 4) 

2 

Where A^ ■ inner jet area * 0.8789^ in 

« fuel flow rate from rotameter readings (CFM) 

5. Outer Jet Mean Axial Velocity (ft/ sec) 

Uo Q^/A^ « 5. 602 ( /HTTP) ^ (A. 5) 

2 

Where A outer jet area » Sir in 
o 

6. Overall Combustor Mean Axial Velocity (ft/sec) 

= .1910(Qj^ + + Qf) (A. 6) 

Where - area of the combustor = 4ir in^ 

7. Inner Jet Equivalence Ratio 

- 0,07501(Q^/m^)/(F/A)g^. (A. 7) 

Where (F/A)^^ is the stoichiometric fuel-air 
ratio “ 0.105 for methane 
“ 0.042 for propane 
p = 0.07501 (Ibm/ft^) for air at STP 

8. Overall Equivalence Ratio 

4>oa “ 0*07501 [Qf/Cm^ + 


f 
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C PROGRAM AVO FIN 
C 

C THIS PROGRAM CALCULATES THE MASS FLUX WEIGHTED AVERAGES 
C FOR THE X<I.J> VALUES FOR AN OPERATING CONDITION 
C X CORRESPONDS TO 03. N2. CO. FUEL. C03. NO. NOX. T AND 
C H20 FOR I"t TO 9 RESPECTIVELY (X(9. J> IS CALCULATED 
C IN THE PROGRAM) J VALUES DENOTE THE RADIAL LOCATIONS 

C AT 05 CM intervals WITH U«1 FDR R«0. 0 AND J«Il FOR 

C R«5 0 CM FOR I«1 TO 5 THE X VALUES ARE SPECIFIED AS 

C A FRACTION OF THE RESPECTIVE MAXIMUM VALUES XMAX. 

C WHICH ARE SPECIFIED IN PERCENTAGE DY VOLUME NO AND 
C NOX ARE SPECIFIED IN PPM AND T IN K 

DIMENSION R< 11 ), U( I I), X<9. 1 1 ) . PROD< 1 1 ) , XMAX ( 5 > . RHOUR < 1 1 > 

1 . SUm<9). AVG(9) 

INTEGER FUEL 

LOGICAL NAMER( 13). NAMEU< 13) . NAMEX (13 ) . NAMEXMU3 > . RUN(2S) 

C 

C ASSIGN LOGICAL UNIT NUMBERS TO FILE NAMES AND READ 

C FROM FILES 

C 

WRITECS. 1 > 

1 FORMAT!' INPUT R VALUES. TYPE FILE NAME R DAT') 

READ <5. lOONAMER 

100 FORMAT ( 13A1 ) 

gall ASSIGNO. NAMER) 

READOi *MR< 1 ). r«l. 11) 

WRITE (5. 2) 

2 FORMAT!' INPUT U VALUES. TYPE FILE NAME C*#«U. DAT') 

READ! 5. lOOINAMEU 

C ALL ASS I GN ! A , NAMEU ) 

READ (A,*) !U!I). I»l. 11) 

WRITE! 5. 4) 

4 F0RMAT!2X, 'INPUT X!IiU> VALUES. TYPE FILE NAME C<>##X.DAT') 
READ!5. 100)NAMEX 

CALL ASSIGN! 1. NAMEX) 

DO 3 I»l. 8 

READ!1. t) !X(I. J). J«l. 11) 

3 CONTINUE 
WRITE!S. 5) 

5 FORMAT!' INPUT XMAX ARRAY. TYPE FILE NAME C*#+XM DAT') 

READ! 5. 100)NAMEXN 

call ASSI0N!2, NAMEXM) 

READ!3. «> !XMAX!I >, 1=1, 5) 

C 

C READ DATA FROM TERMINAL 

C 

WRITE !5, 14) 

14 FORMAT!' INPUT INLET FUEL<PPM> AND FUEL , 1=C3H8, 2=CH4 ' > 

read !5.«) F. FUEL 

WRITE (5, 15) 

15 FORMAT!' INPUT RUN CONDITION') 

READ !5, 151) RUN 



no nr>^ nooono^ m o* ooo>- ►- •-*- nooo— rsooi 


154 


ORIGINAL 
Ol' POOR 


page u 

quality 


iSl 
16 


FORMAT <25A1 ) 

WRITE (6/16) RUN 
F0RMAT<20)(, 2HA1/ /) 


CALCULATION OF MOLE FRACTIONS ON A DRY BASIS 


DO no 1*1,5 
DO no J*l, 1 1 
XU. U)*X(1, U) kXMAXU )/100 
10 CONTINUE 


CALCULATION OF WATER, MOLE FRACTIONS ON A WET BASIS 
AND CONVERSION OF ALL MOLE FRACTIONS TO PPM 


DO 130 U*l. n ' 

IF (FUEL. EQ. 2) GO TO 110 
X<9, J)» 4 /3. #(X(3, J)+X(5, J) ) 

GO TO 119 

le X<9, U)* 2 *<X<3, U) ♦•X<5/ J> > 

19 DO 120 1=1. 5 

X< I, J>= XU, J) #1000000. /( 1+X<9, J) > 

20 CONTINUE 

X<9, J>= X<9, J)#1000000 /U+X(9/J> > 

30 CONTINUE 

OUTPUT OF MOLE FRACTIONS (IN PPM > / VELOC I TY AND TEMPERATURE. 
WRITE(6, 6) 

FORMAT U X, ' R ',2X,-' U ',3X.' 02 ',3X, ' N2 ',2X.' CO ',2X, 

1 'FUEL'/ 2X, ' COO ',2X, ' NO ',1X/' NOX ',1X, ' T 
'/ IX/ 

2 ' H20 /, IX/ 'CM '/ 2X/ ' M/S '/ 3X/ ' PPM '. 3X, ' PPM 

3 OX/ ' PPM '/ 2X/ ' PPM'/ 2X, ' PPM ', 2X, ' PPM '. IX, 

4 ' PPM '/ IX/ ' K '/ IX, ' PPM '. /) 

DO 7 J=l/ 11 

WRITE(6, B)R(J)/U(J),(XU,J),I = 1,9) 

FORMAT ( IX/ F3. 1, 2X/ F5. 2/ IX, 2( IX, F7. 0), IX, F5. 0, 2( IX/ F6. 0)/ 

1 IX, F6. 2/ 2X, F6. 2, 2X, F5. O, IX, F7. 0) 

CONTINUE 

CALCULATION OF AVERAGE VALUES. 


DENOMINATOR 
DO 9 1=1,11 

RHOUR(I) * R(I)#U(I>/X(S, I) 

CONTINUE 

INTEGRATION 

SIMPSON'S RUL.E FDR INTEGRATION BETWEEN R=0. 0 AND R=5. 0 
CALL S I MPS ( RHOUR, SUMO ) 

TRAPEZOIDAL RULE BETWEEN R= 5, 0 AND 5, 1 WITH VALUE AT 5. 1 EQUAL 
TO ZERO 

SUMO = SUMO #0 05»RH0UR( 11 > 
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c. 

C NUMERATOR 

C 

DO 11 Inlf? 

DO 10 J«li 11 
PROD(J)«RHOUR( J> 

10 CONTINUE 

C INTEGRATION 

CAUL S I MRS ( PROD< SUM U > ) 

SUM< I >«SUM( I )+0 05«PR00U1) 

C 

AVO< I >»SUM<nXSUHO 
C 

11 CONTINUE 
C 

C OUTPUT OF AVERAGE VALUES. 

C 

WRITE! 6, 12) <AVG{ I >, 1*1, 9) 

12 F0RMAT(/<6X< 'AVG=',2X,2(1X,F7 0), 1 X, F5 0, 2< 1 X, F6. O) , IX, F6 2, 

1 2X, F6. 2, 2X, F5, 0, IX, F7. 0> 

C 

C CALCULATION OF EFFICIENCIES 

C 

EFF1= (1. -AVG(4)/F)«100. 

C 

IF (FUEL. EG. 2) GO TO 20 
C 

EFF2» ( 1. -(AVG(4)'fAVOO)*67. 63/4B4 22>/F)*100. 

GO TO 30 
C 

20 CONTINUE 

EFF2* ( 1. -'(AVG(4>+AV0(3)*67. 63/189. 76>/F)#100. 

C 

30 CONTINUE 

C 

C OUTPUT OF EFFICIENCV VALUES. 

C 

WRITE (6,17) EFF1,EFF2 

17 FORMAT! /,6X, 'CHEMICAL EFFICIENCY (UNDURNT FUEL ONLY) «',F7. 3, 

1 '7.', /, 6X, 'CHEMICAL EFFICIENCY (UNBURNT FUEL AND CD)«',F7.3, 

2 "/.',//) 

C 

STOP 

END 

SUBROUTINE SIMPS(Y,SUM> 

DIMENSION Y(ll) 

SUM=0 

DO 10 1 = 1/4 

10 SUM = SUM+4*Y ! 2*1 )4-2*Y( 2*1 + 1 > 

H=0. 5 

SUM«H/3«!Y! 1 )+SUM+4«Y! 10)+Y! ID) 

RETURN 

END 
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PROGRAM EFF FTN 

*«##*«»«»*««■»«« 

THIS PROGRAM CALCULATES THE THERMAL EFFICIENCY FOR A 
PROPANE FIRED RUN THE VARIADLE NAMES AND THE DATA 
INPUT FORMAT ARE THE SAME AS FOR AVG. FTN. 


DIMENSION RdD.UUU, X<9, n>, XMAX ( 5 > , PROD ( 1 1 > » H< 1 1 ) 

LOGICAL NAMER< 13). NAMEU< 13), NAMEX( 13>. NAMEXM< 13), RUN( 13) 

ASSIGN LOGICAL UNIT NUMDERS TO FILE NAMES AND READ 
FROM FILES, 

WRITE<5, 1) 

1 FORMAT!' INPUT FILE R. DAT ' ) 

REAP (5, lOO)NAMER 

100 FORMAT! 13A1) 

CALL ASSIGN! 1, NAMER) 

READ! 1, #) !R!I), 1 = 1, 11 ) 

WRITE!5, 2) 

2 FORMAT!' INPUT FILE C»##U. DAT ' ) 

READ! 5, 100)NAMEU 

CALL ASSIGN!2, NAMEU) 

READ!2,#)!U!I), 1 = 1, 11) 

WRITE! 5, 3) 

3 FORMAT!' INPUT FILE e<Ht*X.DAT') 

READ! 5, lOONAMEX 

CALL ASSIGN!3, NAMEX) 

DO 4 1 = 1, 8 

READ!3, #) !X!I, J), J=l, 11) 

4 CONTINUE 
WRITE! 5, 5) 

5 FORMAT!' INPUT FILE C*»*XM. DAT') 

READ! 5, 100)NAMEXM 

CALL ASSIGN!4, NAMEXM) 

READ!4,*)!XMAX!I), 1=1, 5) 

C READ DATA FROM TERMINAL, 

WRITE!5, 6) 

6 FORMAT!' INPUT RUN CONDITION') 

READ! 5, 100)RUN 

WRITE!5, B> 

8 FORMAT!' INPUT ! INLET TEMP )/ 1000, FUEL MASS RATE, P IN MM OF HO') 
READ! 5, *)T1, FM, P 
C 

C COEFFICIENTS FOR CP FOR 02, N2, CO, C3HB, C02, H20 

C 

Al=7. 361 141 
Bl=~5. 3695B9 
Cl =20. 54179 
Dl=-25. 86526 
El = 15. 94566 
Fl=-4. B5B89 
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G1«=0. 5661501 
A2-7. 709920 
D2— 5. 503897 
C2<13, 12136 
D2*~ll. 67^55 
E2*0. 233997 
F2«*~l. 1731Q5 

02- 0, 103B83 
A3-7, 012249 
B3—6. 66B293 
C3-17. 28296 
D3—17, 2B709 
E3-B. 860125 
F3-2. SJ4019 

03- 0, 24477B5 
A4=6. BOOB 
D4-2B. 71 
C4-62. 349 
D4— 107. 19 
E4-69. 802 
F4— 21. 101 

04- 2. 4476 
A5-4. 324933 
D5-20. B0895 
eS— 22. 9459 
D5-16. 844B3 
E5— 7. 935665 
F5-2. 121672 

05- 0. 2408713 
A9-7. 9BB86 
D9— 1, 506271 
C9-6. 661376 
D9— 4. 65597 
E9-1. 696464 
F9— 0, 3706212 
09=0. 03992444 

CALCULATE MOLE FRACTIONS ON A DRY BASIS, 

DO 9 1 = 1,5 
DO 9 J- 1 , 1 1 

X< I, J>=X( I, J)*XMAX(I >/l00, 0 
CONTINUE 

CALCULATE WATER AND MOLE FRACTIONS ON A DRY BASIS. 
DO 10 J=l, 11 

X<9, J>=4. 0#(X(3, J) + X(5, J> >73 
DO 20 I = 1 , 5 

X(I, J)=X( I, J>/< 1+X(9, J) ) 

20 CONTINUE 

X(9, J)-X(9, J)/< 1 + X<9, J) ) 

10 CONTINUE 
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WRim<6. gOOJMUN 

200 FORMATS-'. lOX. 'THERHAU EFFICrENCV', lOX. 'RUN S I3AI///10X, ' R ' 
» I 5Xi ' H S lOX. ' C S3X. ' PROD '/) 

C 

C EVALUATE INTEGRAND AT EACH RADIAL LOCATION 

C 

DO 12 J»li 11 
T-'XIB, Jl/lOOO 0 
C 

C COEFFICIENTS FOR CP VALUE AT THE RADIAL LOCATION J 

A-X< 1» U>»A1 + X(2i U)*Aa+X<3, J)#A3+X<4. U)«A4’^X(5. U)#A3-*-X(9< JMA9 
D«<X(1. J)»D1 + X<£!. J)»D2+X<3. J)»n3+X(4. J)«D4^XOi J>*D0+X(9. J)*D9) 

1 /2 

C»< X(1. J)*Cl-*-X(2iU>*C2-»-X(3. J>*C3+X(4. J)#C4+X(5. J)#C5+X(9i J)*C9) 

1 /3 

D«(X( 1. J)#D1+X(2. JMD2+X(3< J)*^D^+X(4. J>«D4+X(5* J)tD5+X(9. J>#D9) 

1 /4 

E«(X( 1. J)«E1*-X<2, J>#E2-^X(3, J)#E3+X(4. J)*E4+XC&i J)*E5+X<9< J)«E9) 

I /5 

F«(X< 1. U)*F1'^X<2. U)*F2+X«3. J>»F3+X(4. J)*F4-fX(5. J)#F5+X(9. J)»F9) 

1 /6 

0=(X( 1. J)»01-fX<2. JHG2-fX(3i J)^^03+X(4. J)*04-^X(5- J>»C5+X(9i U>*G9> 

1 /7 

C 

C SENSIBLE ENTHALPIES AT THE LOCAL AND INLET TEMPERATURES. 

HT2«T-tHA'fT#{D+T*(G'fT*(D+T»<E'<-T»<F+T*Gn > > ) > 
HTl»TI»(A+THKD-fTI*(C-fTI*<D+Tl*(E+TH^(F+TH»G) ) > > > > 

C 

C increase in SENSIDLE ENTHALPY AT RADIAL LOCATION J. 

H<U)“!(HT2“HT1 > 

C 

C INTEGRAND (CONSTANT FACTORS ARE MULTIPLIED LATER) 

PROD! J)»R( J)-t*U( J)#H(J)/X(e, U) 

WRITE(6, 11 )R(J). H(J>. C, PROD(J) 

11 F0RMAT(10X< F3. 1. 5X. F7 4, lOX. F7 4» 3X< F7 5) 

12 CONTINUE 
C 

C EVALUATE INTEGRAL OVER THE RADIUS 

C 

CALL SIMPS(PROD. SUM) 

C TRAPEZOIDAL RULE FROM R=5. 0 TO R»5. I CM 

SUM=SUM+0. 05«PRDD(ll) 

C 

C MULTIPLY CONSTANT FACTORS. 

C HE« INCREASE IN SENSIBLE ENTHALPY FLUX AT THE EXIT. 

C HE=SUM#100*(P/R)»2*PI«SUM»P»0 01007. R=62400 MM. HG-CC/GM MOLE-K 

C P IN MM. OF HG . 100 CONVERTS M/S TO CM/S. 

C 

HE=0 01007^»SUM*P 
C 

C HEATING VALUE AT INLET TEMP TI 

C HV= VALUE AT 5TP. Q= VALUE AT TI 

C (HV IN KCAL/GM MOLE. Q JN KCAL/GM) 
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HV-484. 22 

A«S»Al+A4-3»A5~4#rt9 
D*<5»Dl+n4“-3«n5-4«H9)/2. 0 
C»=(5*Cl+C4-3#C5-4«C9)/3. 0 
Db( 5*D1+D4-3#D0-4*D9)/4. 0 
E«(S*El+E4-3*<E5-4*E9)/5 0 
Fb<5*FJ+F4“3»F5~4«F9>/& 0 
0»(5*01+04-3#G5-4«G9)/7. 0 

HT2=TI«<A'fTI#<D+TI*<C+TI*{D+TI«-(E+THKF+TI#0> ) ) ) )) 

T«0, 293 

HTl«T*<A+TJ*<DfTI#(C+TI#(D+Tl«(E+THKF+tI»G) > ) > > > 
Q»<HV+HT2-HTl)/44. 0 

CALCULATE THERMAL EFFICIENCY <7.>. 

EFF«lOO<tHE/<Q*FM) 

WRITE<6. I3>SUM( HE, Q 

FQRMAT< 'O', lOX, 'SUM«',F6 4, 5X, 'HE=',F7 4, 5X, 'Q=',F6 3/) 

WRITE<6, 14)EFF, TI,FM, P 

FORMATUOX, 'EFF= ', F7. 3/// lOX, 'TJ*',F5. 3,3X, 'FM»',F6. 4, 'P«',F6. 2> 

STOP 

END 


SUBROUTINE SIMPSCY, SUM> 

C SIMPSONS RULE FOR 1 1 POINTS) STEP SIZE*0. 5 
DIMENSION Y<in 
SUM=0. 0 
DO 10 1 = 1,4 

10 SUM=SUM+4*Y(2*I )+2*Y<2*l+l ) 

H=0. 5 

SUM=H/ 3# ( Y < 1 > +SUM+4# Y < 1 0 ) +Y ( 1 1 > > 

RETURN 

END 
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PROCRAH EFPM. F7N 

««»»«««»«•««#«#» 

THIS PROGRAH CAUCUUATES THE THERHAL EFFICIENCY FOR A 
METHANE FIRED RUN COMMENT CARDS FOR EFF FTN APPLY 
FOR THIS PROGRAM ALSO. 

DIMENSION R(l 1), U( n ). X<9i 1 J )» XMAX < 5) . PROD(l I > , H< 1 1 ) 

LOGICAL NAMER<13>< NAMEU< 13>. NAMEX( 13), NAMEXN( 13). RUN< 13) 
WRITE<55, 1) 

FORMA'-{' INPUT FILE R DAT') 

READ (8, J OO)NAMER 
FORMAT U3A1 > 

CALL ASS10N( li NAMER) 

READ< I, *) <R(1 ), I»l, 11 ) 

WR1TE<5, 2) 

FORMATS' INPUT FILE C#*PU DAT') 

READ (5, 100)NANEU 
CALL ASSIGN(2, NAMEU) 

READ<2, ♦XU(l), I-l, 11) 

WR1TE(5, 3) 

FORMAT!' INPUT FILE C»*»X DAT') 

READ (3, 100)NAMEX 
CALL ASSIGNO, NAMEX) 

DO 4 1 = 1, B 

READO, *) <Xa, U), U=l< 11) 

CONTINUE 
WR1TE(5, 5) 

FORMAT!' INPUT FILE C»#»XM. DAT') 

READ! 5, lOONAMEXM 
CALL ASSIGN!4, NAMEXM) 

READ!4, *) !XMAX!1 >, I«l, !3) 

WRITE!5, 6) 

FORMAT!' INPUT RUN CONDITION') 

READ!5, 100)RUN 
WRITE!5, 8) 

FORMAT!' INPUT ! INLET TEMP )/ 1000. FUEL MASS RATE, P IN MM OF HG ' ) 
READ15. *)TI, FM, P 

COEFFICIENTS FDR CP FOR 02, N2. GO, CH4, COD, H20 
Al«7. 361141 
Bins. 369589 
Cl*20. 54179 
Dl=-25, B6526 
El = 15. 94566 
Fl=~4. 85889 
G1=0 5861 '301 
A2-7 709928 
DD-s-S 503S97 
C2^13 12136 
D2«‘-U , 67955 
E2=5 233997 
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F2«-l. 173105 
02»0 103003 
A3«7. 012247 
D3-"6. 668293 
C3«17. 20296 
03*- 17 20709 
E3>0. 060125 
F3*2 314019 
03*0 2447705 
A4*7 918404 
D4*-ll. 41722 
C4»63. 73457 
D4*~75. 25691 
E4«43. 29269 
F4*-12 56732 
04*1, 49695 
A5*4. 324933 
D5»20. 00095 
C5— 22. 9459 
D5»16. 04403 
E5»-7. 935665 
F5«2. 121672 
05*-0. 2400713 
A9«7. 90806 
D9*-l, 506271 
C9»6. 661376 
D9*~4. 65597 
E9*l, 696464 
F9*-0. 3706212 
09*0 03992444 
DO 9 1=1, 5 
DO 9 0*1, 1 1 

X< I» J>“X< I, J>#XMAX< D/IOO. 0 

9 CONTINUE 

DO 10 J=l, 11 

X<9, U)«2*(X(3, J)+X<5, J) ) 

10 CONTINUE 
WRITE <6, 200) RUN 

200 FORMAT(///, lOX, 'THERMAL EFFICIENCY', lOX, 'RUN: 13A1///10X, ' R ' 

1 , 5X, ' H ', lOX, ' C ', 3X, ' PROD '/> 

DO 12 U=l, 11 
T=X(B, U)/1000. 0 

A=X( 1, J)#Al+X<2, J)#A2+X<3, J)*A3+X<4, J)*A4+X(5, J ) *A5+X ( 9, U ) )^A9 
0*<X(1, U)*B1 + X<2, U)*B2-*-X(3, U)#D3+X<4, J)#I34 + X<5, J>»B5+X<9, U)#D9) 

1 /2 

C*( Xn, U>»C1+X(2, J)*C2+X <3, U)*C3+X (4, J)#C4+X<S, J>*C5+X(9, J)#C9> 

1 /O 

D*<X< 1, J)»D1+X(2, J>*D2+X<3, U)#D3+X(4, J)#D4+X<5, J)*D5+X(9, U)#D9) 

1 /4 

E=<X( 1, U)<fEl+X(2, J!#E2+X<3, U)*E3-4-X<4, J > frE4+X < 5, J ) #E5+X ( 9, J)«E9) 

1 /5 

F*<X< 1, U)#Fl-f-X<2, J)*F2tX<3, J>«F3+X<4, J ) #F4+X ( 5, J ) #F5+X ( 9, U)»F9) 

1 /6 


t 
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0»(X< 1. J)*014.X<2, J)»024‘X<2i J>#034-X(4 i J>#Q4+X(5, J)«09+X(9» J)*09) 
i /7 

HT2«T# ( A+T* ( ntT* < C+T » < D+T* ( < F+T#0 )>)>)) 
HTl-Tl»<A-m»{D+TI#(C+THKDtTI*(E+Tl#(F+Tl#C) ))) ) > 

H(J)*=(HT2-HT1 > 

PROD(J)»-R(J)ftU<vJ)*H<J)/X<e. v>) 

WR1TE(&, inR(vJ>»H(J).C.PRaD(J) 
i 1 FORMAT< 10X» F3. I, 5Xr F7 4, IOXj F7, 4< 3X. F7. S) 

12 CONTINUE 

CALL SIMPS<PRDD, SUM) 

SUM-SUM+0 OS»PRDD< 1 1 ) 

C HE«SUM»100*P/R#2»PI-SUM#P*0 0l007i R«62400 MM, HO-CC/OM MOLE-K 

C P IN MM. OF HO. 

HE=0, 01007»SUM«P 
C HEATING VALUE AT TFiMP T1 

HV-1B9. 76 
A*2»Ai+A4-A5“2«A9 
D»(2*Dl+D4-BS-2*09)/2 0 
C«<2*Cl+C4-CS'-2»C9)/3 0 
D»(2*D1+D4-D5-2»D9)/4 0 
E*(2«El+E4-E5-2*E9)/5 0 
F*{2*Fl+F4-F5-'2#F9)/6. 0 
G=(2*Q1+GA-GS-2#G9)/'7. 0 

HT2*=TI»(A+TI#<D+TI»(C'^TI*(D+Tlt(E'^TI*<F+Tl»0> ) ) ) ) ) 

T»0 293 

HTl“T*(A+TI#(D+Tl*(C+TI»(D+TI*(Ei-THMFTTI*C) > > > > ) 

G=«(HV+HT2-HT1 )/16 0 
EFF"100»HE/<Q#FM) 

WR1TE(6. 13)SUN, HE. Q 

13 FORMAT! 'O'. lOX. 'BUM-'. F6 4, 5X. 'HE-', F7. 4, 5X, 'Q«', F6. □/) 

WRITE<6, 14)EFF. Tl. FM. P 

14 FORMAT! lOX, 'EFF- F7. 3/// lOX. 'TI-'. F5. 3, 3X, 'FM»'.F6. 4. 'P»'.F6, 2) 
STOP 

END 

G 

C 

SUBROUTINE SIMPS!Y, SUM) 

C SIMPSONS RULE FOR U PDlNTSi STEP SIZE-0 5 
DIMENSION Y!ll ) 

SUM«0 0 
DO 10 1-1.4 

10 SUM«SUM+4*Y!2*1 I+E^YIEtl+l ) 

H-0. S 

BUM-H/3» ! Y ! 1 ) +SUM+4*Y ! 10 ) +Y ! 1 1 ) ) 

RETURN 

END 



